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ABSTRACT 


The  observational  and  theoretical  worl- 
on  the  low  frequency  electromagnetic  oscillatior  ^ 
of  the  earth -ionosphere  cavity,  known  as  the 
Schumann  resonances,  is  reviewed.  Additional 
theoretical  considerations  are  introduced  chat 
allow  the  models  to  incorporate  all  the  important 
complications  of  the  real  ionosphere.  The  compari¬ 
sons  between  the  theoretical  predictions  and  the 
observational  data  lead  to  implications  about  Che 
structure  of  the  very  low  ionosphere,  and  the 
geometry  of  ionospheric  perturbations  which  are 
discussed. 


IMTHODUCTlOif 


In  r«o«nt  ymn  a  groirlnc  intATMt  hM  bMn  dAVcIoplng 
in  th«  gAophjrsioal  iaplioAtiona  of  th«  AlAotroignAtic  «*aso* 
nanc'^a  of  th*  cavity  fonMd  batwaan  tha  aax*th  and  ita  iono- 
sphara*  Tha  data  that  can  ba  darivad  from  paaaiva  obaarra* 
tiona  of  thaae  raaonanoaa  uccltad  by  tha  naturally  occuring 
lightning  diachargaa  proaiaaa  to  add  considerably  to  our 
knowladga  of  tha  propartiai  of  tha  icnoaphara.  Thara  are 
savaral  faaturaa  of  thaaa  obaarvations  that  sake  tham  a 
uniqua  way  of  atudying  tha  ionoaphara.  Pirat  of  all>  the 
frequenoiast  which  start  at  about  8  cycles  per  second^  are 
■uch  lower  than  tha  fraquanoias  otherwise  available  for 
asking  ionoepharic  radio  aasauraBants  and  therefore  they  are 
strongly  influenced  by  lower  regions  of  tha  ionosphere  that 
are  difficult  to  study  by  tha  standard  rocket  or  ground 
baaed  tachniquas* 

Secondly^  tha  isaasureaants  involve  world  wide  properties 
of  the  ionosphere  as  tha  waves  that  are  observed  have  trav-^ 
ailed  around  the  earth  aayaral  tiaas,  so  that  all  sections  of 
tha  cavity  can«  in  a  sanseji  be  monitored  from  a  single  obser- 

vll4b^8ita*  A  distribution  of  observing  sites  is  necessary^ 

*  _  . 

how^wr*  to  uniquely  separate  the  geographical  distributions 

f 

of  these  properties*  The  aaasuranonts  can  also  usually  be 
Bade  continuously  so  that  it  is  possAbla  to  investigate  some 
of  tha  ti*ansiant  affects  that  involve  tha  ionosphere* 

Tha  earlier  work  on  this  problem  is  reviewed  in  the 
first  section*  Tha  alactromagnatic  propartima  of  the  iono- 
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sphere  are  examined  briefly  in  the  second  section.  The  large 
varlatione  of  these  properties  with  geographic  position 
necessitates  an  extension  of  the  earlier  theoretical  work 
in  order  to  make  quantitative  use  of  the  data.  A  technique 
to  hrnfiie  this  pi’oblem  is  discussed  in  the  third  section. 

In  tne  1  ourth  section  we  compare  the  results  of  theoretical 
models  v?lth  the  existing  observational  data  and  consider  some 
oi  the  geophj'sical  implications  of  these  comparisons.  Some 
further-  nethematical  details  are  to  be  found  in  the  appendix. 
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'I.  HISTORICAL  REVIEW 

The  first  ideas  that  an  ionosi^ere  was  present  above 
the  earth's  surface  and  capable  of  trapping  electromagnetic 
energy  was  put  forward  in  1902  by  Heavyside  and  Kennelly. 
Viatson  in  two  famous  papers  (Watson,  19l8>  1919)  developed 
'ome  of  the  most  important  mathematical  tools  for  dealing 
with  a  spherical  wave  guide,  and  a  tremendous  amount  of 
literature  is  devoted  to  the  propagation  of  electromagnetic 
waves  through  an  earth-ionosphere  wave  guide  (see  Review 
article  by  Wait,  19^3).  It  is  not  UTitil  1952,  however,  that 
we  find  any  literature  concerning  itself  with  the  resonance 
of  the  entire  wave  guide  system.  Schumann  is  credited 
with  first  studying  the  theoretical  aspects  of  this 
problem  and  also  with  Koni^  of  first  attempting  to  measure 
the  resonant  frequencies  (Schumann,  1952,  195^).  For  this 
reason  the  phenomenon  of  the  earth- ionosphere  cavity  reso¬ 
nances  are  known  as  the  Schumann  resonances. 

The  resonant  frequencies  are  mainly  determined  by  the 
size  of  the  earth  and  the  speed  of  light,  but  the  ionospheric 
properties  modify  these  frequencies  to  a  certain  extent.  The 
lowest  resonant  modes  are  at  sub' audio  frequencies  and  far 
below  the  frequencies  ordinarily  used  for  radio  communica¬ 
tion.  These  frequencies  lie  at  the  low  end  of  the  range 
of  frequencies  that  the  radio  people  designate  as  E,L,p. 


Lightning  discharges  provide  a  broad  spectrum  of  energy  at 
these  frequencies  and  it  is  relatively  easy  to  measure  th-s 
background  energy  at  any  time.  The  spectrum  oi  the  light¬ 
ning  source  is  modified  by  the  resonance  properties  of  the 
earth-ionosphere  cavity  system,  and  this  signature  should 
be  observed  in  the  resulting  noise  signals.  The  number  of 
liglitnlng  strokes  occurring  at  any  one  time  is  in  general 
too  high  to  be  able  to  distinguish  at  these  frequencies 
the  wave  forms  of  individual  flashes,  and  one  must  resort 
to  measuring  the  power  density  spectrum  of  the  background 
noise.  It  was  not  till  Balser  and  Wagner's  measurements 
(Balser  and  Wagner,  i960,  I962,  1963)  that  adequate  analysis 
techniques  were  used  to  extract  the  resonance  Information 
from  this  background  noise.  By  digitizing  two  hours  worth 
of  background  noise  received  on  a  vertical  antenna  and 
numerically  filtering  the  digital  data  they  were  able 
to  Identify  the  first  five  resonant  modes  and  to  estimate 
the  Q  of  these  modes  from  the  shape  of  the  power  density 
spectrum.  These  preliminary  results  are  shown  in  Figure  1, 
Balser  and  Wagner  also  had  attempted  to  determine  some  of 
the  higher  modes  (Balser  and  Wagner  J.  R.  MBS,,  i960)  but 
their  results  were  less  statistically  significant.  They 
have  continued  this  work  using  both  digital  and  analogue 
techniques  and  have  studied  the  time  variations  of  the 
power  in  the  resonant  modes  and  of  the  frequencies  of  these 
modes . 


The  successful  work  of  Balser  and  Wagner  really  marks 
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the  beginning  of  a  great  interest  in  this  phenomenon  which  is 
being  pursued  both  theoretically  and  experimentally.  (Benoit, 
1962;  Chapman  and  Jones,  1964;  Galejs,  I96I,  1962;  Gendrin  et 
Stefant,  1962;  Harris  and  Tanner,  I962;  Pierce,  1963;  PoIk  and 
Pitcher.,  1962;  Raeraer,  I96I;  Row,  I962;  Rycroft,  I963;  Thompson, 
1962;  Wait,  1964.) 


Most  of  the  theoretical  efforts  have  been  aimed  at  im¬ 


proving  the  mathematical  models  so  as  to  better  represent  the 
actual  ionosphere.  Schumann,  in  his  pioneering  work,  showed 


that  the  resonant  frequencies  for  a  perfectly  reflecting  iono¬ 
sphere  would  be  given  by  resulting  in  fre¬ 

quencies  of  about  10.6,  18.4,  26.0,  35*5t  and  41,1  for  the 


first  five  modes.  Schumann  also  corwldered  the  effect  of  a 


finite  conducting  ionosphere  and  showed  that  the  resonant 
frequencies  would  be  lowered  and  the  Q  reduced.  Oalejs  con¬ 
sidered  the  effect  of  an  exponential  conductivity  profile  in 
the  ionosphere  and  showed  that  important  energy  losses  occurred 
in  the  very  low  regions  of  the  ionosphere.  Thompson  considerec’ 
the  influence  of  a  vertical  magnetic  field  and  showed  now  the 
resulting  anisotropic  conductivity  profoundly  modified  the 
cavity  properties,  especially  at  night.  Wait,  Galejs,  and 
Harris  and  Tanner  have  made  contributions  toward  treating  the 
inhomogeneities  of  the  wave  guide  system.  None  of  these  ear¬ 
lier  treatments  have  involved  the  full  complexity  of  the  ac¬ 
tual  ionosphere. 


Other  theoretical  efforts  have  been  aimed  at  studying  the 
effect  of  the  lightning  source  statistics  and  their  geographic 
distribution  on  the  resulting  observed  noise  'oectra  (Raemer, 


noise.  Results  reported  by  Bolser  and  Wagner  I960 


Power  Density  of  Vertical 
Electric  Field 

.fi  o>  on 

m  o  o  o  o  o 
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1961;  Holzer,  1958;  Pierce,  I963).  There  is  some  disagree¬ 
ment  about  the  details  of  the  source  spectra*  Lightning 
discharges  are  complicated  phenomenon  with  several  current 
surges  of  different  time  durations  (Watt,  I96O;  Pierce,  1963)» 
Most  of  the  surges  are  rapid  enough  at  the  frequencies  of 
interest  here  to  be  considered  as  impulse  functions,  but 
the  correlations  between  repeated  strokes  modify  the  result¬ 
ing  spectra  and  enhance  the  low  frequency  end  of  the  spec¬ 
trum.  Also,  some  long  discharges  are  slow  enough  to  limit 
the  energy  in  the  frequencies  of  the  higher  resonant  modes. 
Since  the  observed  spectra  depends  on  the  cavity  properties 
and  the  source  locations  as  well  as  the  source  spectra, 
it  is  not  a  trivial  problem  to  resolve  the  source  spectra 
from  the  observed  Schumann  resonances. 

The  experimental  work  is  simple  in  principle,  as  the 
measurements  are  easy  to  make  with  our  modern  electronic 
technology,  but  a  great  deal  of  analysis  is  required  to 
extract  the  information  available  in  these  meisureraertts , 

The  frequency  analysis  of  the  observed  noise  background  is 
often  done  by  digital  computations,  but  this  is  restricted 
to  spot  data,  as  synoptic  studies  would  be  too  Expensive. 
Analogue  techniques  using  banks  of  filters  must  be  resorted 
to  for  extensive  studies  of  the  Schumann  resonances,  A 
popular  way  of  presenting  such  analogue  results  are  the  time- 
frequency  plots  that  the  sonographs  and  other  similar  instru¬ 
ments  can  produce. 


Observations  are  usually  made  with  vertical  antennas  or 
with  induction  coils.  The  vertical  antennas  are  omnidirec¬ 
tional  receivers,  while  the  induction  coils  can  be  oriented 
to  receive  a  particular  component  of  the  magnetic  field.  The 
vertical  antennas  have  the  advantage,  however,  that  they  only 
respond  to  the  TM  modes  which  is  the  mode  of  the  Schumann 
resonances.  Signals  that  originate  in  the  ionosphere, 
because  of  the  reflection  properties  of  the  air  layer,  come 
down  to  the  earth’s  surface  essentially  in  the  TE  mode. 

These  signals  will  not  participate  in  the  Schumann  resonances, 
and  could,  if*  they  are  large  enough,  mask  the  resonance 
signals  in  the  horizontal  magnetic  field  components.  The 
induction  colls  do  have  the  advantage  that  they  are  more 
easily  shielded  from  the  weather,  as  rain  or  dew  can  short 
out  the  electric  signals  even  when  the  antenna  is  enclosed 
in  a  radar  dome.  (Balser  and  Wagner,  1963).  Because  of  the 
large  vertical  electric  gradient  in  the  air,  there  is  always 
also  the  possibility  that  electric  signals  will  be  produced 
by  modulations  of  the  electrical  conductivity  of  the  air. 

This  does  not  appear  to  mask  the  electromagnetic  signals. 

Even  the  presence  of  local  thunder  activity  does  not  mask 
the  Schumann  resonances  if  the  energy  of  local  lightning 
flashes  is  eliminated  from  the  recordings  by  clipping 
(Wagner-personal  comraunlcatlon) . 

Horizontal  electric  measurements  can  also  be  used,  and 
they  are  closely  related  to  the  horizontal  magnetic  measure¬ 
ments,  They  are  perhaps  the  easiest  measurement  to  make. 
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but  unfortunately,  the  most  susceptible  to  man  msde  Interfer¬ 
ence.  A  very  bad  noise  source  are  the  subhamumlcs  of  the 
power  line  frequency.  These  peculiar  signals  are  apparently 
due  to  non  linear  interactions  between  the  generators  of  the 
power  system,  and  unfortunately  these  frequencies  fall  right 
on  top  of  some  of  the  Schumann  resonances.  They  are  distin¬ 
guished  by  the  fact  that  they  have  a  line  spectra,  but  are 
not  quite  exact  subharmonics ,  and  when  displayed  on  an  oscil¬ 
loscope  the  60  cycle  wave  form  is  observed  to  drift  slowly 
through  the  subharmonic  wave  form. 

The  experimental  programs  are  still  in  the  early  stages, 
but  several  results  iiave  been  reported.  The  original  fre¬ 
quency  structure  of  the  Schumann  resonances  sbsl  reported  by 
Balser  and  Wagner  have  been  substantiated  with  minor  changes 
by  many  observers.  (Qendrin,  1962;  Benoit,  1962;  Chapman, 
1964;  Lokken,  1961;  Rycroft,  1963)  The  genei’al  character  of 
the  amplltute  pattern  has  been  explained  in  terms  of  locali¬ 
zation  of  most  of  the  lightning  activity  in  the  equatorial 
continental  regions  at  local  afternoon  times.  (Balser  and 
Wagner,  1962),  Systematic  shifts  of  the  resonance  frequen¬ 
cies  that  repeat  on  a  dally  basis  have  also  been  reported 
(Balser  and  Wagner,  1962).  These  variations  amount  to  about 
5  %  of:  the  frequency.  Some  frequency  shift  is  to  be  expected 
as  the  source-receiver  geometry  changes  throughout  the  day, 
but  a  systematic  change  in  the  cavity  system  cannot  be  ruled 
out  because  of  the  displacement  of  the  magnetic  pole  from 
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the  rotational  nole.  At  present  this  data  ia  only  reported 
for  New  England  so  that  it  is  difficult  to  untangle  the  two 
effects. 

Sudden  shifts  of  the  Schinann  resonances  amounting  to 
5-7^  also  have  been  reported  inraediately  following  the 
Johnson  Island  high  altitude  nuclear  explosion  of  July  9, 
196^.  These  shifts  were  reported  from  both  New  England 
(Balser  and  Wagner,  I963)  and  Prance  (Qcndrin  and  Stef ant, 
1962}  and  there  is  no  doubt  that  a  systematic  change  in  the 
earth- ionosphere  cavity  was  involved, 

A  summary  of  some  of  these  observations  is  given  in 


Table  I 


TABLE  I 


OBSERVATIONAL  DATA  CQNCERNim 


Mode  number  1 

frequency  in  cps  8.0 

Q  4 

diurnal  variations  of 

f^  in  cpa  to ,2 

transient  shift  of  fj^  in  cps 
following  July  9,  i960 
nuclear  explosion  as 
observed  in  France  -0,5 


SCHUMANN  RESONANCES 

234 
14.0  20.0  26.5 

556 

to.3  to. 5  to. 5 


-0.5 


-0.9 


-0.9 


TABLE  I 


OBSERVATIONAL  DATA  CONCERNII^ 


Mode  number  1 

frequency  in  cps  8,0 

Q  4 

diurnal  variations  of 

ff.  in  cps  .0,2 

transient  shift  of  fj.  in  cps 
following  July  9f  I960 
nuclear  explosion  as 
observed  in  Prance  -0,5 


SCHUMANN  RESONANCES 

2  3^ 

14,0  20.0  26.5 

5  5  6 

to, 3  to. 5  to. 5 


-0.5 


-0.9 


0.9 


31  LOW  FREQUENCY  EI£CTROHAGNETIC  PROPERTIES  OF  I0NC6PIIERE 

The  strong  electromagnetic  influence  of  the  ionosphere 
arises  from  the  presence  of  free  electrons  and  ions.  This 
ionization  results  from  solar  and  cosmic  irradiation  and  is 
the  subject  of  intense  study  these  days,  A  general  under¬ 
standing  of  the  many  factors  that  are  involved  has  developed, 
but  there  is  still  a  lack  of  adequate  quantitative  data, 

(see  Review  article  by  Reid,  1964)  and  many  problems  remain 
to  be  resolved.  We  are  not  primarily  involved  in  this  study 
with  these  problems  of  ionization,  but  eventually,  of  course, 
all  the  indirect  methods  of  probing  the  ionosphere  must  lead 
to  consistent  models  of  this  ionization.  Later  we  shall  see 
that  the  Schumann  resonance  data  imposes  strong  limitations 
on  the  ionization  levels  that  can  be  existing  in  the  lower 
D  or  what  is  sometimes  called  the  C  region.  In  this  section 
we  wish  to  examine  the  behavior  of  sub-audio  frequency 
electromagnetic  waves  in  the  environment  of  the  atmosphere 
and  ionosphere. 

The  electromagnetic  waves  of  the  Schumann  resonances 
have  the  interesting  property  of  exhibiting  large  parts  of 
the  entire  spectrum  of  electromagnetic  wave  beliavior  as  they 
propagate  throug^h  the  atmosphere  and  ionosphere,  starting  as 
classical  electromagnetic  waves  in  the  atmosphere  and  ending 
as  magnetohydrodynamic  waves  in  the  P  regions  and  above. 
Because  of  this  one  must  keep  track  of  many  ionospheric 
parameters,  and  the  classical  subdivisions  of  the  ionosphere  j 
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only  on  the  electron  concentrations  lose  their  significance. 
Strlctp.y  speaking,  the  waves  are  not  electromgnetic  waves 
but  plasoa  waves  and  involve  the  partial  pressures  of  zhe 
plasioa  constituents  as  well  as  the  electric  and  nagnetic 
vectors.  Plasiaa  waves  have  been  intensively  studied  during 
the  past  ten  years,  and  one  can  now  draw  on  a  good  under¬ 
standing  of  their  neneral  properties.  In  the  earth *8  iono¬ 
sphere  the  acoustic  waves  have  velocities  so  slow  compered 
to  the  electromagnetic  waves  that  they  are  largely  uncoupled, 
and  we  can  essentially  ignore  the  pressure  terms  and  consider 
the  waves  as  simply  electromagnetic  waves.  Also,  because 
the  ionosphere  is  a  temperate  plasma  and  our  wavelengths  are 
long  compared  to  the  critical  microscopic  lengths,  the  be¬ 
havior  of  these  waves  is  adequately  described  by  Maxwell’s 
equations,  A  complication: arises  from  the  anistropy  of  the 
ionospheric  electric  conductivity  due  to  the  presence  of  the 
earth's  magnetic  field.  Since^at  the  Schumann  resonance 
frequencies ^conduction  currents  dominate  over  the  displace¬ 
ment  currents  above  about  60  KM,  the  anisotropy  of  the 
electrical  conductivity  profoundly  alters  the  electromag¬ 
netic  wave  behavior  beyond  this  altitude.  The  subject  of 
plasma  electromagnetic  waves  is  extensively  treated  in  many 
elcellent  texts  (Ginzburg,  Allis,  Spitz)  and  numerous 
articles,  and  we  shall  only  briefly  review  the  points  of 
interest  to  our  problems. 

To  solve  for  the  electrical  conductivity  a  relationship 
is  derived  between  J”  and  E  from  the  linc'^rized  equations 
of  motion  for  the  plasma  constituants .  In  these  equations 
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the  frictional  effects  are  usiially  represented  In  terms  of  an 
effective  collision  frequency  between  various  constituents 
iBultiplled  by  their  relative  mean  velocities.  This  concept 
is  not  really  correct  when  the  collision  frequency  depends 
on  the  velocity  of  the  individual  particles  (Ginzburg,  Katz) 
but  the  errors  Involved  are  within  a  factor  of  two  and  can 
usually  be  absorbed  by  making  the  effective  Collision  fre¬ 
quency  depend  on  the  probing  electromagnetic  wave  frequency. 


The  solution  of  these  linear  equations  can  be  simplified 
by  making  use  of  the  large  ion  to  electron  mass  ratio  to 
drop  small  terms.  Depending  on  where  in  the  calculations  these 
approximations  are  made,  different  algebraic  formula  are 
obtained  which  are,  however,  numerically  equivalent.  The 
roost  commonly  used  formulas  are  given  by 

eg  =  - - -  (1) 

^6  (  ) 


Oi  r  Ne' 


Pik  -  i  w 


^  cry  w  — 


(2) 


(3) 


These  equations  assume  that  the  plasma  is  electrically 
neutral  and  contains  no  negative  ions. 


0^  is  the  conductivity  iKirallel  to  the  magnetic  fields 

(Tl  the  conductivity  perpendicular  to  the  field,  and  Ox  the 

Hall  conductivity,  all  given  in  mhos  /neter 

is  the  radio  wave  frequency  in  radians  per  second 

cJe  ia  the  electron  gyro  frequency  in  radians  per  second 

CJi  is  the  ion  gyro  frequency  in  radians  per  second 

N  is  the  electron  and  positive  ion  concentration  In 

Uti  Wen  the  effective  collision 

^  ^  frequencies  between  electrons  and  ions, 

electrons  and  neutrals,  and  Ions  and 
neutrals  respectively  in  stc^ 

A  numerically  equivalent  expression  that  is  more  useful  for 
our  purposes  is  given  as 


0“  _ )iL - 


C^e. 


(4) 


g-  I  +  WJI 


(5) 


From  (4)  and  (5)  we  see  that  tiiree  distinctly  different 
regions  exist  with  which  we  associate  three  different 
regions  of  electromagnetic  wave  behavior. 


(a) 

liZ|  «  1 

(rx«  01  ^  o; 

(b) 

|W| » 1,  (ill«l 

QL«  O'x  «  % 

(c) 

|w)  Jill  »  1 

Ox  02.  Oq 

I  ia 

This  representation  emphasizes  the  role  of  the  collision 
frequencies  in  determining  the  wave  types.  At  the  lower 
Schumann  resonances  CO  is  much  less  than  and  somewhat 
less  than  CJl  ,  so  tliat  we  pass  from  region  a  to  b  to  c  as 

and  Ovn.  become  less  than  and  CJl 

respectively.  (The  ion  gyro  frequency  is  about  30-50  cycles 
per  second  in  the  ionospheric  regions  of  interest  to  us  here.) 
Tables  II  and  III  give  mid- latitude  estimates  of  these 
conductivities  at  7«5  cps  and  include  the  displacement 
contribution. 

In  order  to  study  the  electromagnetic  v/avc'  types  a 
plane  wave  solution  to  Ifexwell's  equations  is  sought.  When 
W  «  1.  the  conductivity  is  isotropic  and  we  have 


^  ijx. 


(6) 


This  expression  is  valid  up  to  heights  of  about  60  KM 
in  the  ionosphere.  At  these  lower  ionospheric  levels  we 
cannot  use  (l)  to  calculate  because  the  ion  popu¬ 

lations  are  much  greater  than  the  electron  populations,  and 
their  contributions  to  Oo  cannot  be  ignored.  In  fact, 
below  about  40  KM  In  the  day  and  50  KM  at  night  the  electri 
cal  conductivity  .is  predominantly  ionic. 

As  long  as  l^coUl  \<^<  the  phase  velocity  of 

the  waves  are  essentially  equal  to  the  speed  of  light  and 


TABLE  II 


HID~IATITUD£  DAYTIME  COMPLEX  CONDUCTIVITY  PROFIUS 

Electrical  Conductivity  at  7.5  cpa  Including 
Displacement  Contribution  in  mKoSjrtneter 


10 

18 

27 

■39 

45 

52 

60 

69 

78 


-4.21  X  10*^0 

(,01  -  4.21)  X  10-10 

(.06  -  4.21)  X  10-10 

-10 


(T, 


(.26  -  4.21)  X  10 

(.43  -  4.21)  X  10 

(1.9  -  4.21)  X  10 

(1.1  -  .041)  X  10 

5.6  X  IQ-^ 

-5 

2,0  X  10 
o  ..  -I  /\**4 


-10 

-10 

-8 


-4.21  X  10-10 
(.01  -  4.21)  X 
(.06  -  4.21)  X 


10-10 
10^10 

(.26  4.21)  X  10’^° 

(.43  -  4.21)  X  lo“^° 

(1.9  -  4.21)  X  lo“^ 

(9.9  -  .421)  X  10“^ 

-7 


2 .0  X  10 

4.8  X  10 


-7 


6.8  X  10-i^ 

7.3  X  10-16 

1.3  X  lo"^^ 

1.3  X  10“^^ 

5.6  X  lo”^^ 

-12 

9.2  X  .10 

2.9  X  lO"^ 

2.7  X  lO”*^ 

3.0  X  10"^ 

-r. 


100 

2.4  X 

10-2 

(1.0  + 

.iu)  X  lo''^ 

1.7  X 

- 

1 

o 

113 

(8.8  + 

.021)  X  lo"^ 

('1.7  + 

.221)  X  10'^ 

2.6  X 

-ii 

10  ■ 

127 

(3.7  + 

.021)  X  10"^ 

(2.1  + 

.141)  X  10’^ 

{3.1  - 

.341) 

142 

(1.8  + 

.021)  X  10* 

(2.8  - 

1.4l)  X  lO"^ 

(.59  - 

1.21) 

18 


TABLE  III 

MID-IATITUDE  NIQHTIMS  COMPLEX  CONDUCTIVITY  PROFILE 

Electrical  Conductivity  at  7,5  cpa  Including 
Displacement  Contribution  in 


in  Kwi 

^0 

'-'J- 

8 

-4.2i  X  10“^° 

-4.21  X  10*10 

2.0 

14 

-4.2: 

'l  3 

c  10“^° 

-4.2i 

.  a 

c  10-“ 

9.1 

21 

(.02 

- 

4.21) 

X 

10-10 

(.02 

- 

4.21) 

X 

10-“ 

4.8 

28 

( .06 

4.21) 

X 

10*10 

(.06 

4.21) 

X 

io‘^° 

3.0 

36 

(.17 

4.21) 

X 

lo'io 

(.17 

4.21) 

X 

10-^^ 

2.4 

45 

(.31 

- 

4.21) 

X 

10“10 

(.31 

- 

4.21) 

X 

10-^° 

1.4 

55 

(1.0 

> 

4.21) 

X 

10*10 

(1.0 

- 

4.21) 

X 

10-10 

3.7 

56 

(^-5 

- 

.421) 

X 

10**^ 

(2.0 

- 

.421) 

X 

10*9 

1.0 

77 

2.5 

X 

lO*'*^ 

(1.5 

- 

.041) 

X 

10*® 

6.0 

90 

3.9 

X 

10-5 

(4.2 

- 

.041) 

X 

10-8 

1.1 

105 

2.4 

X 

10'3 

(7.3 

.161) 

X 

10*^ 

8.4 

120 

1,1 

X 

«»o 

10 

(2.9 

+ 

.231) 

X 

10-8 

(8.0 

137 

(6.1 

+ 

.221) 

X 

10“^ 

(4.5 

- 

1.11) 

X 

10*® 

(^1 

156 

(2.5 

+ 

.351) 

X 

10"^ 

(1.5 

> 

2,71) 

X 

10-8 

(-4.8 

177 

(4.6 

+ 

1.21) 

X 

10“1 

(.68 

tm 

2.91) 

X 

10*® 

(-6.5 

200 

(8.1 

+ 

3.81) 

X 

10*1 

(.28 

- 

2.91) 

X 

10*® 

(-6.7 

226 

(1.1 

+ 

1.4l) 

X 

10* 

(.17 

- 

3.01) 

X 

10"® 

(-6.7 

253 

(.54 

+ 

2.21) 

X 

10* 

(.10 

- 

3.21) 

X 

10*® 

(-6.8 

284 

(5.3 

+ 

8.11) 

X 

10* 

(.03 

- 

1.61) 

X 

10*® 

(-3.2 

0} 

X  10 
X  10 

X  10 
X  10 
X  10 
X  10 


-17 

-17 

-16 

-I5 

-14 

-12 


X  10 
X  10*9 
X  10*® 

X  10*f 
X  10*^ 

-  .191)  X 

-  1.8i)  X 

-  7.11)  X 

-  3.1i)  X 

-  l,2i)  X 

-  .731)  X 

-  .421)  X 

-  .111)  X 


-6 


10 
10"® 

lO**^ 

10-7 


10 

10 

10 

10’ 


-7 

-7 

-7 
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the  wave  fronts  will  be  moving  horizontally.  A  small  amount 
of  damping  due  to  the  ccmduotivity  term  can  bsoome  important, 
however,  becaiise  of  the  large  distances  that  the  waves  are 
travelling.  Above  about  50-60  KM  the  conductivity  term  be¬ 
gins  to  dominate  over  the  displacement  cui*i*ent  term,  and  the 
decreasing  phase  velocity  forces  the  wave  normals  to  become 
more  nearly  vertical.  The  damping  is  much  heavier,  but  the 
waves  travel  only  a  relatively  short  distance  through  this 
highly  damped  region.  As  the  waves  move  slightly  higher  in 
the  ionosphere  the  electron  collision  frequency  decreases 
enough  to  allow  the  electron  motitms  to  become  circular,  and 
the  electrical  conductivity  becomes  anistropio.  It  is  nec¬ 
essary  in  order  to  study  the  wave  behavior  to  develop  the  full 
expression  for  the  propagation  constants. 

llie  expression  for  the  general  case  is  very  complex, 
but  \ie  can  recover  a  simple  expression  making  two  more  approxi¬ 
mations  that  hold  for  most  of  the  Schumann  resonance  waves 
above  about  50  KM.  If  we  ignore  displacement  currents  and 
if  we  make  the  quasi-longitudinal  approximation  0^  Cose}}(^sifie 
where  9  Is  the  angle  the  propagation  vector  makes  with 
the  magnetic  field,  we  obtain  as  a  solution 


(7) 
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Above  about  80  KM  is  so  much  greater  than  (JL  that  the 
quasl-loneltuiiinal  approximation  holds  for  the  Schumann 
resonance  waves  every^^here  except  right  on  the  rnagnetic 
equator.  At  lower  altitudes  the  approximation  is  valid  in  a 
emaller  region  within  30*  latitude  of  the  magnetic  poles. 

Prom  Tables  II  and  III  we  see  that  (Tyc  »  <^1  in  the 
region  between  75  and  110  KM  and  then  we  can  set 


Cose 


(O 


These  are  the  waves  we  usually  associate  with  whistlers. 
The  mode  associated  with  the  top  branch  is  essentially  an 
evanescent  wave.  It  is  often  called  the  ordinary  wave  in 
the  whistler  notation^  and  is  polarized  opposite  to  the 
electron  potation.  The  other  branch  gives  us  the  whistler 
node,  which  tends  to  be  guided  by  the  magnetic  field  and 
propagates  with  a  polarization  in  the  same  sense  as  the 
electron  rotation.  This  wave  propagates  with  little  damp¬ 
ing  in  this  region  of  the  ionosphere. 


Above  l40  KM  we  find  the  ion  collision  frequencies  are 
small  enou^  compared  to  the  ion  gyro  frequencies  to  have  CQ  >0^ 
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Under  these  conditions 


i^co(Jl/^s^e 


(9) 


These  can  be  rewritten  using  (1)  and  (5)  to  give  the 
more  familiar  form  when  iA'n  ^ 


Hiese  are  the  Altei  modes  that  arise  when  the  sound 
velocity  is  assumed  negligible*  The  first  branch  is  often 
called  the  ordinary  wave  in  the  magnetohydrodynamic  notation 
although  it  is  the  same  branch  that  was  called  extraordinary 
in  the  nagneto*ionic  notation.  A  less  confusing  notation  is 
to  call  the  top  branch  the  fast  Alfven  wave>  and  the  bottom 
branch  the  slow  or  guided  Alfven  wave*  Both  waves  propagate 
with  little  damping,  and  the  guided  wave  is  strongly  guided 
by  the  magnetic  field.  At  the  Schumann  resonance  frequencies 
Ox  cannot  really  be  ignored  compered  to  (JL  and 
some  modification  of  the  pure  Alfven  wave  behavior  results , 
These  waves  still  propaiate  with  little  damping,  however, 
beyond  about  I50  KM* 


In  the  region  between  110  and  140  KM  the  wave  types  are 
neither  Alfven  like  nor  whistler  like,  Allis  and  Stix  refer 
to  this  region  as  the  region  of  the  lower  hybrid  resonant 
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frequency.!  but  this  terminology  loses  its  meaning  when  the 
collision  frequencies  are  important  parameters  as  they  are 
in  the  ionosphere  at  tlie  Schumann  frequencies.  This  region 
ropi'esents  another  zone  of  intense  damping.  ®i1b  region  also 
.iccoi’jitr  for  most  of  the  damping  of  the  micropulsation 
enore77  tliat  is  observed  on  the  earth’s  surface, 

Tlius  we  find  the  ionosphere  subdividing  itself  into 
regions  according  to  the  electroerngnetic  wave  behavior  at 
the  Schumann  frequencies.  Ihe  physical  parameters  that 
determine  these  regions  are  the  effective  collision  fre¬ 
quencies  of  the  electrons  and  the  Ions  rather  than  the 
electron  densities,  The  overall  behavior  of  the  cavity 
still  depends  on  all  the  parameters ,  The  theory  for  pre- 
dieting  collision  frequencies  is  somewhat  stronger  than 
the  theory  for  predicting  electron  and  ion  densities  as 
the  important  constants  are  better  known.)  so  that  the 
information  derived  from  the  observations  of  the  Schumann 
resonances  will  still  most  profitably  be  used  to  study  the 
ionization  densities. 
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III.  mTIEHATICAL  FCRitJIATION  CF  THE  SCHU^NN  RESONANCE  PRCBIEM 

Oicj  earlier  reported  theoretical  work  on  the  Schumann 
resonance  problem  has  been  limited  to  treating  ionosj^oric 
models  tlriat  varied  v;ith  altitude  but  not  with  geographic 
position.  This  work  liaa  helped  us  to  gain  a  general  picture 
or  bh3  •rfect  of  the  ionosphere,  but  mb  cannot  really  make 
quantitative  use  of  tiie  Schumann  resonance  data  until  our 
Laeor-tical  models  can  include  the  geographical  variations 
CO  'A ell  ao  the  vertical  variations  of  the  ionosphere.  Sug¬ 
gestions  on  procedures  to  be  used  to  attack  this  problem 
hare  been  made  (Buchsbauir.j  1960;  Galejs>  1962;  Harris,  1962; 
Wart^  1964.1).  In  this  section  v/e  wish  to  discuss  a  practical 
method  of  sebting  up  such  models,  and  also  reviev/  the  more 
classical  formulations, 

?he  Schumann  resonance  problem  can  be  approached 
•jlt'iur  ?.s  a  Hreens  function  pi'oblem  in  which  bho  response 
of  the  cavity  bo  a  variable  frequency  source  is  studied,  or 
:u'i  cliaracterxsti.’  .rrcquency  problem.  The  cliaracterlsblc 
frcquo' Dies  are  thOL..-  .reiuencles  (which  in  real  problems 
are  usually  complex  valued)  that  allow  solutions  of  the 
xiu  o  iuabxoiv:  „  no  source  ternit!.  There  is,  of 

('.out'Ju,  vary  close  co!L;ection  between  these  methods,  '.fhe 
bheoroblcal  results  that  are  discussed  in  the  next  section 
v'rc  cerived  from  a  Qr  len-'  function  approach. 
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The  Greens  funetljr.  yrcblein  in  thcsK.  coordinate  s..  stcus 
where  the  equatioiis  arc  separable  oan  be  Jcjiuileci  o'j  straijtht- 
forward  and  well  cocumerj'ced  ei::e';ifun2t-Oi.  veoliiiiv.ues .  (Morse 
and  Feshbach,  1953;  SonLiierfeld,  19^9)  I’ne  procedure  consists 
of  defining  a  surface  of  the  cooralnate  svstein  or  w-iich  the 
source  is  located  aiid  developing  the  source  S4.n.",ularLty  i' 
terras  of  the  eigenfunctions  associatea  i;ith  that  surface.  The 
function  defining  the  behavior  of  the  solutior.  in  the  re¬ 
maining  coordinate  is  not  an  eigenfunction,  oi:t  since  there 
are  always  independent  solutions,  the  roiiiaicin.;  boundary 
conditions  can  be  satisfied  by  a  proper  balance  of  these 
Independent  solutions.  In  many  practical  problems  thei'O 
are  several  ways  of  deflriing  the  surface  so  that  different 
formulations  of  the  problem  result,  and  these  uifferent 
formulations  lead  to  different  physical  interpretations  as 
well  as  different  numerical  procedures.  The  eigenfunctions 
act  like  standing  waves  but  otherwise  they  have  no  wave 
propagation  properties.  In  fact  the  same  eigenfunctions 
often  appear  in  potential  problems  as  do  in  viave  propaj^a  oic- 
problems.  The  propagation  wave  behavior  is  only  seen  lx.  the 
coordinate  direction  that  is  perpendicular  to  the  surface  on 
v;hich  the  source  was  defined . 

In  spherical  problems  the  defining  surfaces  most  often 
used  are  the  surfaces  /?=  fo  and  0=0  (assuming  the 
source  Is  at  the  pole).  In  the  former  v.ase  the  eigenfunctions 
are  S1V\^  COS  (cose)  .  The  vxave  propagation 

behavior  is  exlilbited  in  the  R  direction  by  the  spherical 
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bessel  functions  li^(kR)  arid  ,  In  the  latter  case 

the  eigenfunctions  are  Sin^  Cos 
where 

proper  boundary  behavior.  (Soraerfeld,  1949)  The  wave  propa¬ 
gation  behavior  is  exhibited  in  the  0  direction  by  the 
Legendre  functions  (-Cos 6)^  * 


and  h  i  (kf^)  h  (kR)  , 
are  complex  values  that  give  the  radial  functions 


The  former  approach  using  spherical  surface  harmonics  is 
more  familiar  and  therefore  more  often  used.  However,  in  a 
problem  involving  propagation  of  waves  around  the  earth  it 
would  appear  more  logical  to  use  the  latter  approach,  U'atson 
(Watson,  1919)  showed  how  one  could  transform  a  solution  in 
the  0  harmonics  into  one  in  the  radial  harmonics  ar*d  that 
the  latter  representation  was  indeed  much  more  efficient  for 
many  computational  purposes,  (see  also  Bremmer,  1949)*  The 
radial  eigenfunctions  that  arise  in  this  form  are  the  v/ave 
guide  modes . 


Both  types  of  solutions  may  be  exhibited  Li  a  single 
integral  which  has  the  form  (Wait,  1962;  Bremmer,  1949) 


u 


(-Cos©  )  ill 


(10) 


.-I  _ 

The  contour  C  encloses  the  poles  of  5in  lill  and 
gives  the  spherical  harmonic  series  representation.  The 
contour  can  also  be  deformed  so  that  the  value  of  the  integral 
Is  given  instead  by  the  residues  of  the  poles  of  ^ 
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which  are  the  radial  eigenfunctions  or  the  wave  guide  mode 
solutions.  This  is  the  transformation  that  is  called  the 
V/atson  transformation. 

The  characteristic  frequency  determination  essentially 
merges  these  two  approaches  since  the  characteristic  solutions 
are  eigenfunctions  in  all  the  coordinate  directions ,  (See 
also  Bremmer,  19^9;  V/aitj  1964.)  If  we  consider  the  integrand 
(10)  as  the  frequency  60  is  varied  a  pole  of  at  certain 

.  -t  ^ 

values  of  60  will  merge  with  a  pole  of  Sin  pTT  ,  ie.  at 
x)  =.  n.  .  The  set  of  these  frequencies  are  the  charac¬ 
teristics  frequencies  or  free  oscillations  of  the  system 
and  would  cause  the  integral  (10)  to  diverge  since  a  source 
term  is  implied  in  the  representation.  The  real  part  of 
these  characteristic  frequencies  are  the  resonant  frequencies 
of  the  system,  and  the  imaginary  pert  of  the  characteristic 
frequency  can  be  used  to  estimate  the  Q  of  the  system 


The  meaning  of  the  9  of  a  system  tends  to  be  a  matter 
of  definition.  The  more  rigorous  definitions  involving  the 
rate  of  energy  dissipation  are  not  often  amenable  to  direct 
observation  especially  in  continually  excited  systems,  so  that 
other  observable  quantities  are  used.  'Ihe  most  often  used 
quantity  is  the  width  of  the  observed  resonance  peak  in  the 
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power  specti’uitie  For  simple  oacillating  systems  setting  ^ 

where  Ap  is  the  separation  of  the  half  power  points,  is  con- 

I  I  J  flv,  r 

sistant  with  the  definitions  7^  =  ~  as  well  as 

p.  ^  Jicucle 


9-- 


VFhen  the  Q  is  high  these  definitions  are  also  consis- 
tant  for  any  system.  Practical  difficulties  arise  when  the 
Q  is  low  enough  to  merge  neighboring  peaks  in  the  power 
spectrum.  In  such  cases  the  half  power  widths  give  an  under¬ 
estimate  of  the  other  two  definitions  of  Q  .It  is  very 
likely  that  this  occurs  in  analysing  the  Schumann  resonance 
peaks . 


These  analytic  approaches  that  have  been  outlined  are 
limited  to  models  that  allow  our  equations  to  separate, 

Unf oi’tur.ately,  some  of  the  most  interesting  aspects  of  the 
actual  earth-ionosphere  cavity  are  complications  that  des¬ 
troy  this  siraiplicity ,  The  influence  of  the  magnetic  field 
which  is  so  very  important  in  the  night  time  ionosphere  pre¬ 
vents  the  equations  from  separating  in  spherical  coordinates 
unless  the  magnetic  field  is  vertical,  (Thom|B3ona  I962 , ) 

The  day-night  and  latitude  asymmetries  introduce  further 
complications.  If,  however,  we  are  to  make  quantitative  use 
of  the  observations  of  the  Schumann  resonances  we  must  find 
a  practical  way  of  incorporating  all  these  complications  into 
our  mthematical  models ,  A  brute  force  attack  on  the  problem 
Is  difficult  uecause  of  the  strong  three  dimensional  charac- 
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ter  of  the  problem.  It  turns  out,  hov.'ever,  that  the  problem 
can  be  broken  down  into  two  steps,  each  of  which  is  fairly 
straightforward  and  conveniently  organized  for  act’jal  compu¬ 
tations,  The  technique  is  essentially  based  on  the  use  of 
transmission  line  equivalences  to  waveguide  propagation,  a 
concept  that  is  used  widely  and  very  successfully  by  the 
microwave  circuit  people  in  handling  very  complicated  wave¬ 
guide  problems  (Montgomery,  et  al,  19^8).  The  power  of  the 
method  derives  from  the  fact  that  any  waveguide  mode  can  be 
represented  as  propagation  along  a  transmission  line,  and 
this  representation  is  ejtact  as  long  as  only  one  waveguide 
mode  is  involved.  In  such  a  case  the  waveguide  behavior  can 
be  reduced  to  a  problem  in  circuit  theoi^.  When  other  modes 
are  involved  and  the  modes  couple  together  at  boundaries  or 
interfaces  of  the  system,  some  further  complications  are 
introduced.  The  beautiful  advantage  of  these  methods  in 
the  Schumann  resonance  problem  comes  about  because  at  the 
frequencies  of  interest  only  the  principle  or  TEM  mode  can 
propagate  without  very  severe  damping. 

If  we  consider  the  field  due  to  a  point  source  in  a 
waveguide,  it  can  be  represented  as  a  contour  Integral  in 
which  the  poles  of  the  integrand  represent  the  waveguide 
mode  contributions  and  the  branch  cut  represents  the  essen¬ 
tially  spherically  spreading  wave  energy  and  the  near  field. 
This  breakup  is  often  described  as  the  discrete  and  continuous 
spectra  (Brekhovskiki,  I960}.  The  continuous  spectra  becomes 
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negligible  compared  to  the  waveguide  inodes  at  distances 
large  compared  to  the  waveguide  height  (Pekeris,  19^8). 

When  the  height  of  the  guide  is  less  than  a  half  wave  length 
only  the  zeroth  order  modes  can  propagate  without  severe  damp¬ 
ing.  Por  the  Schumann  resonances  the  wavelengths  are  frac¬ 
tions  of  the  earth's  circumference,  while  the  height  of  the 
guide  is  only  50  to  60  KH,  so  that  the  higher  order  modes  are 
damped  out  in  distances  of  a  few  hundred  kilometers.  The 
energy  that  is  observed  at  the  resonant  frequencies  has 
travelled  around  the  earth  several  times  so  that  the  con¬ 
tinuous  spectra  is  also  negligible.  Thus  the  wave  energy  is 
essentially  confined  to  a  single  mode,  the  TSM  or  zeroth 
order  TM  mode. 

For  each  such  mode  the  relationship  between  the  trans¬ 
verse  components  of  E  and  H  as  described  by  I'fexwell's  equa¬ 
tions  are  of  the  same  form  as  the  relationship  between 
voltage  and  current  in  a  transmission  line.  By  propel’  normal¬ 
ization  of  the  transmission  line  parameters  it  is  also 
possible  to  make  the  energy  flow  along  the  transmission  line 
equivalent  to  the  energy  flow  along  the  waveguide.  If  only 
one  mode  is  involved  a  properly  normalized  transmission  line 
is  an  exact  equivalent  of  the  actual  waveguide. 

Consider  propagation  along  a  guide  in  the  X  direction  of 
a  TEM  mode  with  no  dependence.  Prom  Maxv/ell's  equations 
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we  write 


=  f(r-i££j)Ez 


n  I  ^Jjy 

'  '(c-ietj') 


Along  a  tranamisBlon  line  we  have 


AV 

Ix 


=  -£I 


(12) 


(13) 


(14) 


(15) 


(16) 


^  is  the  series  impedances  and  y  the  shunt 

admittance  per  unit  length.  Waves  propagate  along  such  a 
line  with  a  propagation  constant  given  as 


k= 


(17) 
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The  characteristic  impedance  of  the  line  is  given  as 


z 


cHara.^ 


(17) 


In  order  to  equate  the  energy  flow*  we  set 

.Ml 

V  =  )  Ep 


(19) 


1  =  -H;. 


(20) 


Alt  is  taken  as  the  height  whvre 
disappears  so  that 


essentially 


o  o 

(  Hy  is  essentially  constant  in  the  range 
for  the  T.EM  mode) 


(21) 

Z  =  0,  All 


The  energy  flow  along  the  ground  is  so  small  we  take  the 

bottom  as  on  the  earth's  surface.  In  a  microwave  guide  Alt  is 

simply  the  height  of  the  guide,  but  in  the  ionospheric  guide 

the  upper  walls  are  diffuse.  Because  of  the  exponential 

increase  of  conductivity  with  height,  and  the  accompanying 

rapid  disappearance  of  it  will  be  accurate  enough  to  take 

Alt  as  the  height  where  the  conductivity, ^6CJ  and  assume 

m 

j  cKf  = -t6£0- ALT.  £2.(0) 

o 


that 
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Now  from  (12)  and  (14)  we  can  write 


/A  LI 


(22) 


iV 

i-Y- 


u).  Art  - 


it  a  correction  due  to 


(23) 

fii- 

L  (T-teu  h  ?*■ 


Now  we  can  identify  Y  ^  with  the  parametera  of  the 

waveguide  ae 


(24) 


(25) 


K;(  ia  the  propagation  conatant  of  the  waveguide  BK>de  and 
froai  (17)  *ee  it  will  alao  be  the  propagation  conatant  of 
the  tranaaiaaion  line,  the  correction  factor  !■  taken 
account  of  by  the  value  of  kx  •  a  loaadeaa  guide  the 
TKM  mode  propagatea  with  the  velocity  of  light  and  1?^=  . 

For  such  a  guide  2  repreaented  by  inductive  elementa  and 
y  ia  represented  by  capacitive  elementa.  In  the  actual 
earth  ionoa];^ere  guide «  the  leakage  of  energy  out  of  the 
guide  leada  to  an  added  reaiative  component  in  Z  ,  and  the 
dlsaipatlon  cT  energy  in  the  guide  leads  to  an  added  conduc¬ 
tive  component  in  Y  . 
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The  continuity  conditions  for  V  and  I  arc  the  same  as 
those  for  t  ^  and  ,  so  that  as  long  as  the  waveguide  has 
constant  dimensions  the  analogy  between  the  t^o  systems  .;ill 
be  correct.  Even  when  the  diraensions  of  the  waveguide  change, 
the  analogy  holds  as  long  as  other  modes  are  not  generated  by 
these  changes.  When  sudden  discontinuities  arise  in  the 
guide  it  is  usually  necessary  to  introduce  other  waveguide 
modes  in  order  to  satisfy  the  boundary  conditions.  When  these 
modes  are  non-propagating  their  effect  can  be  represented  by 
an  added  shunt  reactance,  but  to  determine  the  value  of  the 
reactance  one  must  solve  a  boundary  value  problem  (Montgomery, 
19^8),  For  gradual  changes  of  the  waveguide,  this  effect 
is  unimportant. 

The  parameters  of  the  equivalent  transmission  line  also 
depend  on  the  height  of  the  guide,  so  that  changes  in  the 
ionospheric  height  will  cause  Impedance  changes  in  the  equiva¬ 
lent  transmission  line  and  these  changes  v/111  modify  the 
resulting  propagation.  In  order  to  demonstrate  chat  these 
Impedance  changes  correctly  simulate  the  effects  of  an 
actual  gradual  change  of  ionospheric  height,  let  us  consider 
a  tapered  guide  or  a  sectorial  horn  for  which  the  exact 
solutions  of  the  TEM  mode  are  kno’.n  (  Slater,  1959)  •  The 
walls  of  such  a  waveguide  are  formed  by  the  surfaces  d-  6^ 
and  (9  =  in  cylindrical  coordinates.  The  TEM  moae  is 

represented  by 

AT,(ke)^-  BN„(ke) 


(26) 
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L, 


(27) 


Because  of  the  change  of  height  as  we  move  out  along  the 
guide  (height  is  defined  here  as  the  arc  length  along 
P  *  const)  the  equivalent  transmission  line  parameters 
are  given  as 


z  =  2^ 

Y  =  Ya/eA© 


The  transmission  line  equations  are  thus  given  as 


(32) 
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Y»  ^  ^  we  have 

T=  AT„(kf)+  BWo(ke)  (34) 

In  this  case  we  get  the  exact  solution  v/ith  our  trans¬ 
mission  line  equivalence.  This  is  so  because  no  other  wave¬ 
guide  modes  were  needed  to  satisfy  the  boundary  conditions. 
For  a  varying  taper  this  would  not  be  true,  but  as  long  as 
the  change  of  the  taper  is  gradual  the  errors  involved  in 
neglecting  the  higher  order  modes  are  very  small  (Slater, 

1959.) 

Harris  and  Tanner  (Harris,  I962)  proposed  an  interesting 
analytic  technique  to  handle  waveguide  Inhomogenieties,  but 
their  method  does  not  take  account  of  changes  Of  ionospheric 
height  and  thus  does  not  properly  keep  track  of  the  wave¬ 
guide  impedances , 

The  generation  of  other  modes  depends  on  the  mismatch  of 
the  E  and  H  components  at  the  point  where  the  guide  pro¬ 
perties  change,  and  these  mismatches  have  a  (l-CoScTo^  like 
iependance  where  <5*0  is  the  change  of  angle  involved.  Thus 
the  criteria  for  ignoring  mode  coupling  involves  the  wave¬ 
guide  changing  its  snape  slowly  compared  to  the  height  of  the 
guide.  It  is  important  to  note  that  the  transmission  line 
equl*  alence  does  not  depend  on  the  assumption  that  the  wave¬ 
guide  properties  change  slowly  compared  to  a  wavelength. 


36 


which  is  a  necessary  condition  for  the  W.K.B.  approximations 

which  have  been  used  to  treat  V.L.F.  propagation  in  non- 

uniform  waveguides  (Wait,  1964a).  The  wavelengths  of  the 

Schumann  resonance  frequencies  are  actually  much  longer  than 

the  lengths  of  the  transition  zones  separating  the  day  and 

night  ionospheres  for  example,  and  these  changes  car^  therefore, 

set  up  partial  reflections  which  are  accounted  for  by  the 

are 

transmission  line  equivalences,  but^ unaccounted  for  in  the 
W.K.B.  approximations. 

This  all  indicates  that  we  can  accurately  describe  the 
propagation  properties  of  a  section  of  the  earth-ionosphere 
waveguide  for  the  Schumann  resonance  problem  in  terms  of  two 
parameters  such  as  the  characteristic  impedance  of  the  guide 
and  its  propa^tion  constant.  To  determine  these  parameters 
for  a  section  of  the  earth  ionosphere  waveguide  ue  must 
concern  ourselves  with  the  ionospheric  structure  and  the 
magnetic  field  strength  and  its  orientation.  This,  however, 
will  only  Involve  us  with  a  one-dimensional  problem  as  the 
transmission  line  equivalence  allows  us  to  separate  out  this 
section  of  the  waveguide  from  the  rest  of  the  system.  In 
order  to  solve  the  one -dimensional  problem  we  have  recourse  to 
many  techniques .  The  method  used  in  this  study  involved  a 
cylindrical  geometry,  and  the  air  and  ionosphere  were  des¬ 
cribed  by  a  large  number  of  layers  of  constant  properties. 

We  shall  leave  the  description  of  the  method  actually  used 
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to  obtain  solutions  for  the  appendix.  Many  solutions  must 
be  carried  out  in  order  to  get  an  accurate  parameterizatio: 
for  all  the  various  sections  of  the  earth-ionosphere  system. 

We  have  been  using  the  term  “transmission  line"  which 
is  a  one-dlraensional  concept,  but  in  this  particular  problem 
we  are  really  dealing  with  a  two-dimensional  transmission 
surface.  This  means  that  the  surface  impedance  2  must  be 
a  2nd  order  tensor,  the  principle  axes  of  the  tensor  being 
oriented  by  the  earth's  magnetic  field.  Actually,  the  one- 
dimensional  solutions  showed  that  some  TE  mode  vjas  produced 
by  the  twisting  effect  of  the  magnetic  field  on  the  electro¬ 
magnetic  waves  In  the  ionosphere  so  that  2  should  really 
be  considered  as  a  4th  order  tensor.  The  amount  of  energy 
involved  in  this  other  mode  was  very  small,  however,  and 
even  the  anisotropy  of  ^  was  minor  so  that  2  can 
for  most  purposes  be  considered  a  scalar.  The  anisotrop-j  of 
^  is  not  a  complication  as  long  as  our  axes  are  aligned 
with  the  principal  axes  of  2  ^ 

Once  the  appropriate  properties  of  the  transmission  sur¬ 
face  are  determined,  there  still  remains  the  problem  of  solv¬ 
ing  for  the  response  of  this  transmission  surface  to  lightning 
excitation.  This  is  now  only  a  two-dimensional  problem  and  is 
far  simpler  to  deal  with  than  the  original  three-dimensional 
problem.  The  computations  in  this  study  were  done  by  solving 
a  two-dimensional  network  ttot  was  the  lumped  cii.  '.uit  repre¬ 
sentation  of  the.  transmission  surface  equations.  A  schematic 
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representation  of  such  a  network  is  shown  in  figure  2.  The 
methods  of  computation  are  discussed  in  the  appendix. 

It  would  also  be  possible  to  obtain  analytic  solutions 
of  the  transmission  surface  equations  by  perturbation  tech¬ 
niques  j  although  we  did  not  use  this  approach  in  oir  compu¬ 
tations  . 

In  order  to  gain  some  insight  on  the  behavior  of  an 
equivalent  circuit  representation,  let  us  consider  a  uniform 
earth-ionosphere  model.  The  network  representing  the  trans¬ 
mission  surface  will  be  set  up  along  latitude  and  longitude 
lines.  The  impedances  connecting  the  nodes  together  will  be 
proportional  to  the  distance  between  the  nodes  and  inversely 
proportional  to  the  width  of  surface  it  is  representing. 

The  admittances  connecting  the  nodes  to  ground  will  be  pro¬ 
portional  to  the  area  of  surface  associated  with  the  nodes. 


■-  Sin6  4/Ae  (36) 

Y  =  Yo 


(37) 
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If  the  eource  is  at  a  pole  so  that  we  can  ignore  the 
(j)  dependence,  the  network  equations,  in  the  limit  of  an 
infinite  number  of  nodes,  lead  directly  to  the  Legendre 
differ<mtial  eqiuition. 


iV  -  -  Y  R^sik  S  X 
le  1“ 
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41  =  -  4V/sikeA<fi 
IB  / 
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Let  J  =  COS  ©  and  since  ^  -  -  fe* 
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The  resonance  condition  for  this  network  will  occur  for 

1,  2,  3  or  feoK  which  arc  the  correct 

resonance  conditions  for  a  uniform  spherical  waveguide.  The 
changing  impedances  of  the  network  as  we  approach  the  poles 
because  of  the  decreasing  area  are  responsible  for  setting 
up  partial  reflections  that  modify  the  waves  and  cause  a 
phase  stagnation  around  the  pole  producing  a  standing  wave 
pattern.  This  is  one  way  of  looking  at  the  behavior  of  a 
caustic,  but  one  must  remember  that  this  same  phase  stagna¬ 
tion  will  occur  if  the  source  is  placed  on  the  equator,  even 
though  now  at  the  caustic,  which  is  on  the  equator  at  the 
opposite  aide  from  the  source,  the  impedance  values  are  not 
changing  with  position. 

To  test  the  accuracy  of  a  lumped-circuit  approximation, 
the  resonance  frequencies  for  a  circuit  with  nodes  every  10 
degrees  of  latitude  were  compared  with  the  resonance  fre¬ 
quencies  of  a  continuous  system.  The  values  obtained  are 
shown  in  Table  IV.  In  order  to  make  the  calculations  practi¬ 
cal,  a  small  amount  of  damping  was  used. 

The  first  three  resonant  modes  are  accurate  to  within 
or  better  in  the  lumped- circuit  approximation.  Somewhat 
greater  errors  can  be  expected  when  the  source  is  placed  on 
the  equator  and  a  two-dimensional  lumped  circuit  is  used. 
Using  a  grid  of  nodes  every  10  degrees  of  latitude  and  20 
degrees  of  longitude  to  represent  a  homogeneous  ionospheric 
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TABLE  rV 


Resonance  Frequencies  of  Spherical  Waveguide  (almost  lossless) 


Lumped  Circuit  Approximation 

True  Value 

Mode  Number 

NoJes  every  10  degrees 

10.38 

10.55 

1 

18.18 

18.31 

2 

25.6 

25.9 

3 

32.6 

33.3 

4 

39.5 

40.9 

5 

46.2 

48.4 

6 
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cavity  the  frequencies  of  the  first  three  resonant  modes  were 
too  low  by  1,  2.  and  35^.  A  discussion  of  the  error  sources 
is  made  in  the  appendix  but  it  is  safe  to  say  that  the  errors 
will  be  about  the  same  when  the  inhomogeneous  models  are 
studied.  When  these  circuits  are  used  to  study  perturbations 
of  the  cavity  one  can  expect  the  changes  produced  in  the 
resonant  frequencies  to  have  the  same  relative  accuracy. 
Increasing  the  number  of  nodes  in  order  to  increase  the 
accuracy  greatly  increases  the  computational  time,  but  the 
first  order  correction  obtained  by  simply  increasing  the 
frequencies  of  the  first  three  modes  by  1,  2,  and  35^  should 
give  ample  accuracy  for  this  problem  with  the  10*  x  20*  grid 
spacing. 
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IV.  GEOPHYSICAL  IMPLICATIONS  OP  THE  SCHUMANN  RESONANCE  DATA 

The  use  of  the  Schumann  resonance  data  to  study  phenom¬ 
ena.  in  the  ionosphere  is  still  in  the  early  stages^  but  a 
few  positive  results  are  already  beginning  to  accrue.  In  this 
section  we  wish  to  examine  some  of  the  observational  data  in 
the  light  of  the  model  calculations  that  were  outlined  in  the 
previous  section^  and  discuss  some  of  the  implications  that 
arise. 

The  particular  points  to  be  examined  are  the  effect  of 
the  earth’s  magnetic  field,  limits  on  the  lower  D  or  C  layer 
ionization,  amplitude  and  frequency  variations  of  the  reso¬ 
nant  frequencies,  and  transient  effects  of  particular  iono¬ 
spheric  perturbations. 

Magnetic  Field  Effects 

No  data  at  present  is  available  that  shows  the  height 
dependance  of  the  Schumann  resonance  signals,  although  this 
could  be  a  method  of  directly  measuring  the  conductivity  of 
the  lower  D  or  C  regions.  It  is  worthwhile,  nevertheless, 
to  examine  the  theoretical  solutions  as  a  function  of  height, 
in  order  to  gain  physical  insight  on  the  wave  behavior.  The 
solutions  were  carried  out  for  a  layered  model  of  the  atmos¬ 
phere  and  ionosphere,  and  the  electromagnetic  fields  at  each 
layer  boundary  were  computed  as  well  as  the  fields  on  the 
earth’s  surface.  Prom  these  results  the  energy  flows  could 
be  calculated,  and  the  Q's  computed  from  the  rate  of  energy 
dissipation.  In  figures  3  and  4  we  show  the  energy  diBSi|iation 
plotted  as  a  function  of  height  for  the  raid-latitude  daytime 
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and  nlghtlme  model.  These  results  follow  closely  the  expected 
behavior  discussed  in  the  second  section  and  i^raphasize  the 
Important  role  played  by  the  earth's  magnetic  field.  Th» 
three  electromagnetic  windwsj  a  radio  window  in  the  atmos¬ 
phere  below  the  ionosphere,  a  whistler  window  in  the  upper  D 
and  lower  E  region,  and  an  Alfven  window  in  the  P  region  and 
above,  are  clearly  seen.  These  last  two  windows  are  due  to 
the  magnetic  field.  In  the  daytime  model  little  energy 
actually  reaches  the  P  region,  but  at  night  one  third  of 
the  energy  completely  escapes  the  cavity.  This  results  in 
lower  Q^s  for  the  nightime  propagation  even  though  the 
phase  velocities  may  be  higher.  Calculations  have  been 
carried  out  to  show  that  if  any  of  this  energy  is  guided 
back  to  the  cavity,  it  produces  a  negligible  effect  on  the 
resonance  properties. 

The  magnitudes  of  the  electric  and  magnetic  vectors 
show  different  behaviors.  The  radial  electric  field  is 
essentially  inversly  proportional  to  CTjrr  -L6cJ  and 

thus  rapidly  disappears  above  the  effective  ionospheric 
height  which  is  hiuse  taken  to  be  where  (Iq  =  6  CO  .  The 
horizontal  magnetic  wave  field,  however,  remains  quite  con¬ 
stant  through  this  region.  The  field  begins  to  twist  through 
the  whistler  mode  region  and  then  damps  in  the  transition 
region.  The  magnetic  field  that  leaks  through  this  region  in 
the  nightime  ionosphere  suffers  very  little  more  damping. 


Ill" 
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Energy  Dissipation  end  E.M.  Wave  Types  at '7.5  cps 
Day-Time  Ionosphere  Mid  -  Latitudes 


Figure  3  iSnergy  dissipation  o.f  Schumann  resonances 
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The  importance  cf  the  rnagnetic  field  is  further  seen  on 
comparing  the  solutions  derived  for  different  latitude  en¬ 
vironments,  In  these  computations  the  ionization  levels  v;ere 
the  same  for  the  different  models,  excepting  the  polar 
model,  but  the  strength  and  inclination  of  the  magiietic  field 
simuiated  the  actual  field  of  the  earth,  3oma  increase  of 
the  C  layer  ionization  which  is  believed  due  to  cosmic  ray 
ionization  v/as  allov/ed  in  the  polar  model£  ,  The  propagation 
characteristics  of  these  models  are  given  in  ’Table  V,  The 
lAiase  velocities  were  determined  from  the  phase  relation¬ 
ships  of  the  model  and  are  given  as  a  fraction  of  the  velo¬ 
city  of  light,  “ttie  Q.’s  were  determined  from  the  ratio  of 
the  energy  stored  to  the  energy  dissipated. 

The  results  of  Table  V  indicate  that  the  magnetic  field 
variations  with  latitude  produce  changes  in  the  propagation 
characteristics  that  are  as  great  as  the  day-night  variations. 
The  day  and  night  models  have  different  effective  heights, 
however,  which  will  also  modify  the  waveguide  impedances. 

Some  differences  of  the  propagation  characteristics 
were  found  depending  on  the  direction  of  the  propagation, 
but  these  were  not  very  important  and  Table  V  represents  an 
average  of  the  N-S  and  E-W  propagation  characteristics. 

Ionospheric  Profiles 

The  variations  apparent  in  Table  V  indicate  that  only 
crude  interpretations  of  the  resonance  data  can  be  attempted 
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until  the  full  three-dii^Bional  aodelfi  are  worked  out. 
NevertheleBB,  some  modlficationB  of  proposed  iwiospherlc 
models  can  be  made  at  this  stage  because  of  the  large  dis¬ 
crepancies  Involved.  The  Sohuiaann  resonance  data  can  Most 
usefully  be  used  to  study  the  very  low  D  regions  or  the  C 
region  because  at  these  low  frequencl*3  the  conductivity 
becomes  an  Important  parameter  in  the  wave  propagation  at 
much  lower  altitudes  than  is  true  for  V.L.F.  frequencies. 
Also,  since  the  rocket  data  is  not  very  dependable  at  these 
low  latitudes,  and  the  inpoz^tant  ionization  and  recombina¬ 
tion  reaction  rate  coefficients  are  poorly  known,  the  proper¬ 
ties  of  this  region  are  less  well  determined  than  those  of 
higher  regions  of  the  ionosphere.  An  example  of  this  state 
of  affairs  is  shown  in  figure  5  which  gives  daytime  electron 
concentration  estimates.  These  estimates  include  both 
theoretical  determinations  based  on  estimates  of  the  ioniza- 
ulon  xates  and  the  rates  of  the  modifying  reactions,  and 
experimental  determinations  using  ground  based  and  rocket 
techniques.  Most  of  the  recent  theoretical  work  leads  to 
much  higher  ionization  levels  than  the  old  ground  based 
determinations  as  given  by  Pejer  and  Gardner  (Pejer,  1955J  '  ^ 
Gardner,  1934).  The  more  recent  experimental  data  as 
summarized  by  Knetch  appears  to  support  these  higher 
estimates,  (Knetch,  19^3 •)  The  theoretical  work  of  Moler 
and  of  Arnold  and  Pierce  are  exceptions  to  this  trend, 
however,  and  their  estimates  are  much  more  in  agreement  with 
the  older  experimental  values,  (Mol»3r,  I960;  Arnold,  1963.) 
Much  less  data,  but  equally  severe  discrepancies, are  to  be 
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found  in  the  ion  concentration  estimateB,  (Bourdeau,  1959; 
Wlipple,  196^.) 

It  is  clear  from  the  results  ehoiwi  in  figures  3  and  4 
that  the  ionization  levels  at  around  40-50  KM  are  very 
important  for  the  Schumann  resonances,  and  yet  we  find  th»b 
estimates  of  ionization  at  these  levels  varying  by  a  factor 
of  50:1. 

The  Schumann  resonance  data  must,  therefore,  be  able 
to  discriminate  between  the  two  trends  shown  in  figure  5. 

In  figures  6  and  7  are  shown  the  conductivity  values,  G'q  , 
that  result  from  these  profiles.  These  estimates  also 
include  the  ionic  contribution.  The  EWKB  model  incorporates 
Whipple's  theoretical  estimates  of  electron  concentration, 
Knethh's  summary  data  on  electron  concentration  and  Bourdoau's 
measurements  of  ionic  conductivity.  The  same  sources  are  used 
for  the  nightirae  model  NWKB,  The  Dll  model  used  Moler's 
estimates  of  electron  concentration  and  ionic  conductivities 
that  are  a  compromise  between  Bourdoau's  measurements  and  Warn 
measurements  reported  by  Whipple.  (Bourdeau,  1959;  Whipple, 
1964.)  The  S.I.D.  model  is  based  on  estimates  reported  by 
Nicolet  and  Aiken  (Nlcolet,  i960)  and  by  Belrose  and  Cetlner 
(Belrose,  I962).  The  nigiittlme  model  N12  also  relies  on  Moler's 
estimates,  and  uses  the  same  low  level  ion  concentrations  as 
Dll.  The  collision  frequency  values  used  at  these  levels  are 
essentially  those  given  by  lalgarno  (Dalgarno,  I96I) .  This 


Conductivity  Models  (cTo)  Night-Time 


Figure  7  Nightime  conductivity  models  and  resonance  properties 
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data  is  aunnarlzed  in  l^ble  VI. 

The  resonance  properties  that  are  also  given  In  figures 
5  and  7  are  based  on  the  Ionospheric  zaodels  being -uniformly 
distributed  over  all  the  earth's  8urface>  and  are  therefore 
not  realistic,  but  still  the  results  clearly  discriminate 
against  the  DWKB  and  NWKB  models.  The  phase  velocities  and 
Q's  of  these  models  are  so  low  conypared  to  the  data  obtained 
from  the  Schumann  resonances  that  one  must  reject  these 
models  as  providing  too  much  Ionization.  Whipple  discusses 
the  possibility  that  his  estimates  are  too  high  and  suggests 
that  dust  particles  are  responsible  for  sweeping  up  much  of 
this  low  level  Ionization.  It  seems  unlikely  that  the  dust 
levels  are  high  enough  to  accomplish  this  much  of  a  change# 
and  If  so  the  Schumann  resonance  data  should  be  very  sensl*- 
tive  to  Influxes  of  dust  such  as  result  from  the  meteor 
showers,  for  which  there  seems  to  be  no  evidence.  It  Is 
more  likely  that  the  reaction  rate  constants  that  are  Impor¬ 
tant  in  the  theory  tend  towards  the  values  used  by  Moler. 

The  more  complete  calculations  using  the  Dll  and  N12 
models  for  the  day  and  nightlme  ionization  levels  gave  very 
reasonable  results  and  most  of  the  calculations  to  be  dis¬ 
cussed  in  the  rest  of  this  section  used  these  models. 
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In  order  to  teot  the  senBltivity  of  the  Schunum  reeo- 
nancee  to  changes  of  these  ionization  levels «  the  resonance 
properties  of  slightly  perlurbed  aodels  were  also  coinputed. 
These  results  are  given  in  Table  VII.  Ionization  increases  . 
in  the  C  and  lower  D  regions  always  reduce  the  resonant 
frequencies  and  the  Q's«  but  the  effect  is  less  pronounced 
at  the  higher  frequencies.  Ionization  increases  at  higher 
levels  produce  more  conplex  effects.  In  the  day  time  in¬ 
creases  of  the  £  layer  ionization  increases  the  frequency 
and  Q  of  the  first  laode^  but  decreases  the  frequency  and  Q 
of  the  third  mode.  At  night  the  Q's  were  increased  for  all  . 

.  of  the  first  three  modes  when  the  £  layer  ionization  in- 
creasedj  but  the  Q  of  the  first  BK>de  wcs  decreased  and  the 
Q  of  the  third  ande  increased  when  the  F  layer  ionization 
increased.  It  would  be  very  difficult  to  interpret  the 
E  and  P  layer  variations  from  observations  of  the  Schumann 
resonances,  but  there  are  many  other  techniques  for  obser¬ 
ving  these  layers.  The  Schumann  resonances  are  most  sensi¬ 
tive  to  changes  of  the  D  and  C  regions  and  the  use  of  the 
resonance  data  will  probably  be  confined  to  studying  these 
regions  of  the  ionosphere. 

Diurnal  Variations 

The  diurnal  variations  of  the  backgroxmd  E.L.F.  noise 
power  are  perhaps  the  best  documented  features  of  the  Schumann 
resonance  phenomenon.  (Holzer,  1956,  1956;  Aarons,  1956; 

Raeraer,  I96I;  Balser,  1962;  Polk,  I962; 
Stefant,  I963.)  The  first  investigators  realized  that  the 
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general  activity  wee  related  to  the  world  wide  thunder  atom 
activity*  The  ratio  of  the  avplitudea  of  the  various  nodes 
is  nainly  determined  by  the  location  of  the  receiver  relative 
to  the  nodes  of  the  standing  wave  pattema  set  up  by  the 
source*  This  is  modified  by  the  source  spectrum  and  by  the 
complications  of  an  inhomogeneous  wave  guides  and  by  the 
frequency  dependence  of  the  damping  factors*  Raemer  attempted 
to  adjust  the  source  distributions  in  order  to  explain  the 
observed  spectra,  but  because  of  still  unresolved  factors^ 

I 

concerning  the  wave  propagation  complications  and  the  source 
spectra,  these  results  are  soBMWhat  qualitative*  Nevertheless 
the  general  agreement  with  the  known  thunderstom  areas  is 
good  enough  to  warrant  some  confidence  in  our  general 
understanding  of  the  sources  of  the  observed  resonances* 
Similar,  but  more  qualitative,  results  have  been  reported 
by  other  investigators  (Holzer,  1958;  ijalser,  1962;  Rycroft, 
1963 j  Stefant,  1963),  During  the  Northern  Hemisphere  winter 
the  main  areas  of  thunderstorm  activity  occur  on  the 
equatorial  land  areas  at  around  16:00  local  time.  During 
the  summer  important  sources  also  develop  in  the  mld- 
latitude  continental  regions  of  the  Northern  Hemisphere 
(liandbook  of  Geophysics,  i960).  Figure  8  shows  the  diurnal 
power  variations  that  would  be  predicted  for  a  localized 
source  occurring  at  16:00  local  time  on  the  equator  recorded 
in  New  England.  Some  recent  winter  time  data  for  New 
England  provided  by  Balser  and  V/agner  is  also  shown.  The 
general  features  of  the  observed  data  could  be  explained  by 
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E.S.T. 


Diurnal  Power  Variations  1®^  Mode 

Figure  8  Theoretical  and  observed  diurnal  power  variations  of 
first  resonant  mode  for  New  England 
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inox*eAslng  the  loportance  of  the  Afrlomn  aoureee  and  Inolud* 
Ing  a  .background  of  more  aoattered  aourcea.  In  theae  aolu- 

V 

tiona  it  ahould  be  pointed  out  that  not  only  are  the  varla-' 
tlona  of  the  vaveguida  propagation  oonatant  brought  Into  play, 
but  alao  the  iapedanoe  changea  of  the  waveguide*  A  stoat 
ifliportant  paraaeter  in  deteradning  the  inpedance  la  the 
effective  ionoapheric  height^  Alt*  The  ionospheric  height 
plays  a  further  role  when  the  aa^litudes  of  the  vertical 
electric  field  are  being  considered^  since  in  the  trans- 
Bission  line  analogy  is  equated  not  to  V  ,  but  to  V/Alt* 

Added  inforaation  can  be  provided  by  exasiining  the 
power  in  the  various  components  of  the  mgnetic  variations 
produced  by  the  Schuaann  reaonancea*  Unfortunately >  some 
the  reported  aagnetio  data  has  not  specified  the  components 
involved  * 

In  order  to  aake  more  quantitative  uae  of  the  power 
level  data  one  should  really  have  data  from  a  network  of 
stations  distributed  over  the  entire  earth* 

While  there  is  a  good  general  understanding  of  the 
reasons  for  the  observed  diurnal  r^wer  variations,  little 
work  had  been  done  to  explain  the  observed  diurnal  variations 
of  the  apparent  resonant  frequencies*  Systematic  studies 
of  these  variations  have  only  been  reported  by  Balser  and 
Wagner  (Balser,  1962),  but  some  examples  of  these  variations 
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are  also  given  by  Stafant  (Stefant«  1963 )•  These  observers 
have  all  noted  that  the  variations  were  not  sys  tens  tic  anong 
the  various  nodes  and  therefore  were  not  prinarlly  due  to 
systeaatic  variations  of  the  ionosphere.  Balser  has  sug*> 
gested  that  the  variations  are  due  to  the  interplay  of  the 
various  hansonlc  nodes  as  the  source  receiver  geometry  is 
changed.  As  was  discussed  in  the  last  section,  the  wave 
field  in  a  spherically  syneU'ic  waveguide  can  be  described 
in  tarns  of  a  set  of  spherical  hambnics,  (COSO) 

The  magnitude  of  each  term  in  the  harmonic  series  depends 
on  how  close  the  excitation  frequwioy  is  to  the  resonant 
frequency  for  that  harmonic,  and  the  ft.  When  the  Q  is  high 
and  the  excitation  frequency  is  equal  to  a  resonant  frequency 
the  amplitude  of  that  harmonic  is  so  large  the  other  terms 
can  be  ignored.  In  such  a  case  one  would  observe* a  peak  in 
the  power  spectrum  at  this  resonant  frequency  for  a  broad 
frequency  source  almost  everywhere  except  at  the  nodes  of  the 
harmonic  funcclon.  As  the  Q  becones  lower  the  other  terms 
in  the  harmonic  series  become  important  further  away  from  the 
node  zones,  and  also  the  shape  of  the  BOiu*ce  frequency  spectrum 
can  begin  to  alter  the  observed  peak  frequencies.  Also  as  as- 
symetries  are  brought  into  the  waveguide  system,  the  harmcmio 
series  must  be  expanded,  aadh  term  splitting  into  2n  4-  1  terms, 
(Co5B)-^  P^Ccose)^  .  when  the  as\jmmetry  is 

small  the  frequencies  associated  with  the  m  values  for  a 
given  n  value  are  close  together,  and  the  effect  is  known 
as  line  splitting.  In  such  a  case  even  higher  values  of 
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are  neceoBary  In  order  for  the  reaonant  frequency  to 
excite  only  one  hamonlc  node.  With  the  laodest  Q  values 
that  we  find  in  the  Schuzaann  resonances  the  resonant 
frequencies  oust  excite  quite  a  bit  of  energy  in  aany  other 

hsraonic  modes  and  the  line  9plittin&  though  present^  will 

\ 

be  masked.  As  the  receiver  moves  around  rilative  to  the 
source  the  ratio  of  the  amplitudes  of  the  excited  modes 
change  and  these  changes  modify  the  observed  spectrum. 

Another  way  of  describing  this  situation  is  shown  in 
figure  9*  Here  we  show  the  power  levels  for  two  frequencies 
near  the  resonant  frequency.  It  is  noticed  that  the  node 
zones  move  as  the  frequency  varies.  This  would  not  be 
possible  if  only  one  harmonic  mode  was  involved.  It  is 
readily  apparent  that  the  ratio  of  the  power  of  these  two 
frequencies  varies  with  position^  so  that  the  power  spectrum 
shape  will  also  vary.  The  resulting  varlatioiis  of  the 
apparent  resonant  frequency  are  shown  in  figure  10  for  two 
different  Q  values.  Although  the  computations  in  figure  10 
v/ere  done  for  a  uniform  waveguide  model  the  changes  in  the 
apparent  resonant  frequency  for  the  Q  of  5  are  comparable  in 
magnitude  to  the  observed  diurnal  variations.  When  the  more 
realistic  Dll  '■  H12  model  is  used  the  predicted  apparent 
resonant  frequencies  and  their  diurnal  variations  are  in  close 
agreement  with  the  observations  as  is  shown  in  figure  11. 

Again  the  agreement  v/ould  be  improved  if  the  localized 
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source  was  augmented  by  a  baclcground  of  more  widely  scattered 
soiirces.  It  appears^  however,  that  the  essential  features  of 
the  diurnal  variations  are  explained  by  the  mode  mixing,  and 
no  particular  systematic  variations  of  the  Ionosphere  are 
needed.  In  agreement  with  Balser's  suggestion. 

Since  the  shape  of  the  power  spectx*um  Is  used  to  compute 
an  apparent  Q  as  well  as  the  apparent  resonant  frequency,  the 
observed  Q's  should  also  show  a  diurnal  variation.  A  sum¬ 
mary  of  the  predicted  results  from  the  Dll  -  N12  model  Is 
given  in  Table  Vlll,  The  two  dlnnnslonal  model  used  to  cal¬ 
culate  these  results  used  the  propagation  characteristics 
given  In  Table  V,  as  well  as  the  appropriate  effective 
Ionospheric  heights  and  Interpolated  those  values  as  a 
function  of  latitude,  longitude,  and  frequency.  The  longi¬ 
tude  dependence  was  assumed  constant  except  for  a  20  degree 
strip  connecting  the  day  and  night  Ionospheres.  The  Southern 
and  Northern  hemispheres  were  assumed  to  be  symettlc  so  that 
the  tilt  of  the  earth’s  magnetic  axis  was  Ignored,  as  well 
as  any  north-south  assymetrles  produced  by  the  tilt  of  the 
earth's  axis  of  rotation  from  the  ecliptic. 

Transient  Variations 

Perhaps  some  of  the  most  interesting  applications  of 
the  Schumann  resonance  studies  are  their  use  in  tracking 
large  scale  ionosiAieric  perturbations  associated  with  special 
events.  The  only  well  documented  example  la  the  effect  of 
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th«  Starf iBh  high  altitud*  nuclaar  exploalon  of  July  9,  1962 
(Balaar,  I9632  Qendrln*  1962),  It  appaara  that  a  Byateaatlc 
ohanga'  of  tha  raaonanca  propartlaa  oecurrad  iBMdiately  f  ol> 
lowing  tha  axploaion  ovar  Johnaton  Xaland,  and  tha  raaonant 
fraquanelaa  of  tha  firat  four  aodae  war#  daoraaaad  by  •33, 

•33«  *9  and  .9  cpa  (Stafant,  1963),  aa  obaarvad  in  Franca, 
lha  affaet  daeayad  away  in  a  pariod  of  hours.  A  aiiallar 
dacraaaa  of  tha  first  aoda  was  obaarvad  in  New  England 
(Balaer,  1963)*  The  tiaa  raaolution  of  the  effect  la  of 
course  not  sharp  since  it  takas  several  minutes  to  get  a 
statiatically  significant  determination  of  the  resonant  fre> 
quenoies.  Since  the  resonance  properties  represent  a  world 
wide  averagir^  of  the  waveguide  paraaaters,  it  Is  clear  that 
these  anomalies  must  have  bean  causad  by  a  very  extensive 
perturbation  of  the  ionosphex^a.  It  now  appears,  from  radio 
propagation  studies  made  in  Japan  and  in  the  United  States, 
that  the  entire  Bicific  area  was  iwediately  heavily  ionized 
to  low  altitudes  by  the  explosion.  (Nseda,  1964;  Zmuda,  1963) 
These  results,  as  we  shall  see,  are  substantiated  by  the 
Schumann  resonance  data.  This  ionization  occurred  within 
the  first  second  of  the  explosion.  The  time  resolution  of 
the  maasureawaritB  do  not  allow  much  closer  definition,  and  It 
is  still  not  clear  what  were  the  primary  sources  of  this 
ionization.  Some  of  the  Initial  1 >nization  may  have  been 
fast  gsaaa  irradiation.  If  this  were  the  case  the  sources 
of  the  gawm  radiation  must  have  been  shot  up  far  above  the 
explosion  site  since  radio  and  optical  effects  were  observed 
in  the  western  part  of  the  United  States  (Odencrantz,  1962; 
Zmuda,  19631  Holer,  1964). 
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Other  sourceB  of  ionization  were  probably  supplied  by 
the  neutron  decay  products.  Studies  of  the  geometry  of  this 
radiation  indicate  that  the  entire  Iticif ic  would  be  involved 
(Kenney,  1963;  Oelbersann,  1964).  Further  ionization  resulted 
from  the  dumping  of  er  rgetic  particles  that  were  put  into 
trapped  orbits  by  the  explosion  (Zmuda,  1963).*  A  more 
localized,  but  longer  lasting,  source  of  ionization  would  be 
the  decay  of  fission  products  (Williams,  1962). 

The  Schumann  data  is  not  suited  to  tracking  the  rapid 
sequence  of  events  that  followed  the  explosion,  but  it 
does  throw  some  light  on  the  geometry  of  the  ionization. 

In  order  to  study  thin  problem  we  first  considered  as  a 
perturbation  model  heavy  ionization  in  an  area  of  about 
5000  km  by  8000  km,  centered  over  JohnilDn  Island  in  the  night 
side.  The  ionization  model  used  was  similar  to  the  ioniza¬ 
tion  levels  that  occur  in  the  artic  regions  during  PCA  events. 
If  such  ionization  occurred  uniformly  ovvsr  the  entire  earth 
the  first  mode  would  be  lowered  to  3*7  ops,  but  when  the  zone 
was  restricted  to  46  x  10  square.,  kilometers  around  Johnston 
Island  the  lowering  of  the  first  mode  resonant  frequency  was 
only  .03  cps.  In  figure  12  we  show  the  phase  changes  that 
this  perturbation  would  cause  at  7*5  cps  for  a  source  locat¬ 
ed  just  south  of  India.  There  is  some  distortion  of  the 
I^^'se  contours,  but  over  most  of  the  earth's  surface  the 
phases  are  hardly  changed.  Such  a  map,  of  course,  could  not 
be  drawn  up  from  direct  observations  of  the  Schumann  reso¬ 
nances,  since  the  actual  sources  are  random. 


After  Electron  Dumping 
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It  is  evidsnt  from  thsss  rssults  ttmt  ths  ionization 
effacts  of  tha  Starfish  axplooion  wara  far  wra  axtansiva  In 
araa  than  tha  nodal  usad  in  fl^ura  12.  As  a  sacond  nodal  va 
considarad  a  lass  intansa  ionization  aqulvalant  to  tha  S.I.D. 
nodal  of  figure  3$  but  spread  out  over  tha  entire  laeif ic 
araa  and  covering  about  one  third  of  tha  earth* s  surface. 

Tha  phase  changes  of  this  nodal  are  shown  in  figure  13.  This 
nodal  does  produce  effects  large  enough  to  be  conparad  to 
the  actual  affects  of  July  9*  It  should  be  noted  in  figure 
13  that  the  caustic  which  is  narked  by  tha  large  araa  of 
almost  constant  phase  disappears  because  of  tha  savara 
assynetry  introduced  by  tha  boab  perturbation. 

In  order  to  conpara  the  effect  of  this  nodal  with  tha 
observations  nada  in  Prance  and  in  New  England  wa  oust  also 
consider  the  source  locations,  since  tha  apparent  resonant 
frequencies  vary  with  the  latter.  Therefore,  in  figures  14 
and  13  we  plot  the  resonant  frequency  shifts  of  tha  first 
node  that  would  be  observed  in  New  England  and  in  France  as 
a  function  of  the  source  location.  If  the  major  lightning 
sources  M  the  time  were  located  in  Australia  the  observations 
in  France  would  be  predicted,  but  the  observations  in  New 
England  would  be  not  quite  predicted.  Additional  sources 
located  in  central  and  northern  North  America  might  have 
improved  the  fit,  but  it  certainly  seems  that  the  ionization 
model  used  is  neither  too  extensive  nor  too  severe.  The 
areas  that  are  narked  null  zones  represent  areas  that  are 
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nodac  relative  to  the  receiving  site.  One  cannot  really 
carry  this  stiKly  such  further  without  a  better  knowledge  of 
the  source  locations,  and  a  wider  distribution  of  observing 

sites* 


It  is  also  possible  to  study  the  effects  that  are  to  be 
expected  frca  soae  of  the  naturally  occurring  ionospheric 
perturbations.  In  figure  16  we  plot  the  predApted-effect  of 
an  S.I.D.  The  ionization  aodel  was  the  sane  one  used  to 
■odel  the  Starfish  effect,  but  the  perturbatiCH'.  occtirred  on 
the  daylighz  side  not  the  night  side.  In  figure  16  the 
apparent  resonant  frequency  shift  ia  plotted  as  a  function 
of  the  receiver  location  for  a  source  at  16 i 00  local  tiee, 
which  is  depicted  as  in  eastern  South  America.  The  general 
level  of  the  effect  is  comparable  to  the  Starfish  effect. 

In  figure  17  we  show  the  predicted  effect  of  a  P.C.A. 
event*  The  ionization  perturbation  used  is  more  intense 
than  the  3*1.0*  model,  but,  since  the  effect  is  confined  to 
the  auroral  zones,  a  much  smaller  area  is  involved  and 
the  Schumann  resonances  ere  less  effected*  Not  only  are 
the  magnitudes  of  the,  two  effects  different,  but  also  the 
pattern  of  the  effects  are  quite  dissimilar.  It  appears 
that  one  might  be  able  to  interpret  the  geometry  of  the 
ionospheric  perturbations  froai  a  knowledge  of  the  pattern 
oT  the  Schumann  resonance  perturbations.  At  present 
there  does  not  seem  to  exist  enough  data  to  undertake  such 
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a  The  only  recorded  study  ^  such  a  naturally  occur* 

ring  perturbation  was  by  Chapaan  and  Jones  (Chapoan^  19^) 
who  reported  a  substantial  Increase  of  the  first  aode  frequency 
associated  with  a  P.C.A.  event.  This  would  seea  soaewhat 
unlikely j  but  warrants  further  study. 

Additional  inforaation  can  be  obtained  fron  the  aapli* 
tude  of  the  resonance  signals.  Unlike  the  frequency  values, 
the  amplitude  values  can  be  influenced  by  local  perturbs- 
tlons  directly  over  the  receiver  or  the  source.  This  can  be 
a  disadvantage  for  some  studies,  but  an  advantage  for 
other  studies.  An  interesting  application  of  this  local 
effect  has  been  made  in  order  t#  study  the  D  layer  adjust- 
Bents  to  a  solar  eclipse  (Bgmett,  19d{|).  In  such  a  case  the 
vertical  electric  field  is  strongly  inversly  dependant  on. the 
effective  ionospheric  height. 
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SUIMIRY 

In  reviewing  the  past  work  in  this  field  and  in  our  own 
contributions  to  the  subject  we  have  tried  to  show  that  the 
phenonenon  of  the  Schumann  resonances  are  reasonably  well 
understood*  This  does  not  hsply  that  the  subject  is  closed, 
but  rather,  as  we  examine  the  parameters  that  are  closely 
linked  to  the  resonance  behavior,  it  is  evident  that  the 
subject  should  continue  to  provide  important  information 
concerning  properties  of  the  ionosphere  that  are  difficult 
to  Obtain  by  other  means*  The  two  areas  that  are  best 
suited  to  study  by  means  of  the  Schumann  resonances  are 
the  very  low  regions  of  the  D  layer,  and  the  large  scale 
picture  of  the  geographic  distribution  of  the  ionospheric 
properties*  In  order  to  make  fuller  use  of  the  information 
inherent  in  these  signals,  it  will  be  necessary  to  have  a 
broader  geographical  coverage  of  observing  sites  than 
exist  at  present*  With  a  growing  interest  in  this  field, 
it  is  hoped  that  this  situation  will  change. 
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In  ■•etion  III  it  vns  •hown  that  the  full  traataent  of 
the  earth*ionoeph«re  eaTe«t^ide  could  be  Mpneted  into  two 
parte*  nie  firet  pert  ooneiited  of  finding  the  TEN  node 
propegetion  cherecterietiee  of  verioue  segnente  of  the  wave- 
guides  end  the  eeoond  pert  eonsieted  of  fitting  theee  eeg- 
aeftte  together  end  etudying  the  properties  of  the  complete 
spherical  waveguide.  We  wish  to  briefly  outline  here  the 
aethod  used  for  atteoking  the  first  part^  and  In  appendix 
B  we  #11  outline  the  aethod  used  to  do  the  second  part. 

In  a  given  seoticti  tht  ionospheric  parameters  may  be 
considered  as  being  functions  of  altitude  only.  Because  of 
the  eagnetic  fields  hcaeverj  the  equations  still  do  not 
separate  in  spherical  coordinates.  Cylindrical  coordinates 
however^  allow  the  equations  to  separate  and  also  keep  the 
topological  feature  of  the  spherical  surface^ that  the  arc 
lengths  vary  in  proportion  to  R.  Thus  assuming  no  z  depend- 
ance  snd  an  £  dependence  on  0  we  can  write  Maxwell  *8 
equations  as 

iJje  +  fe  -  y?  Er  =  ^2 

Ar  r  r  ^ 


A-3 
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=  -  %6^B  -  %r^r  - 

^7  ^  +^ir^r +^16^6  *■' 

=  (V  fr  +  (Trt  ^6  ^-  t?  ■ 

The  dlaplaceflMnt  contribution^  -UcO  is  included  in  the 
value  of  the  diagonol  eonduetivlty  teroB  / 

Eliminating  £r  and  Hr  froo  theaa  equations  we  can  write  the 
four  remaining  equations  as 


iX 

17 


A-7 


X 


Ee 

H? 

Ei 


(  0  )  (o  )  i  0  ] 

-/(T  (- 

I  (Trr  /  '  r^rr  /  '  oVr 


A-^ 


A-9 


as 


A  convenient  aethod  of  Integrating  theee  equations  involves 
the  use  of  oatrixants  (Fraser,  Duncan  and  Collar,  193B| 
QantaacVtM^  i960).  A  foraal  solution  of  A»7  is  g^iven  as 

Y(r)  -  X(ft.)  A-ii 


If  x(r«)i8  known  or  if  a  coiridination  of  eleaents  of  x(r) 
and  x(r«)  are  known,  a  knowledge  of  allows  one  to  solve 
for  the  renaining  term.  In  a  waveguide  problem  the  boundary 
impedances  give  one  ratios  of  the  eleaents  of  X,  and  then 
A-10  becomes  a  characteristic  value  problea  to  be  solved  for 
the  characteristic  wave  ntambers. 


X.  is  called  the  aatrizant  and  it  nay  be  deteraindd 
by  an  iterative  integral  procedure 


J2^  ^  I  +  I  +  A(n)  A (G)  J-fi  -f 


h 


A-11 


Also  we  have 


If  sources  produce  a  discontinuity  in  X  of  I't  r^  we  have 


Y(r)--  S, 


A-13 
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In  a  region  where  A  is  constant  the  laatrliant  takes  a  simple 
font  which  can  be  verified  by  substitution 


Jtr,  =  6  =  1 4-(n-ra)A  +(tH')  A-A  +  /\-A  Af- 
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A-l4  is  much  more  convenient  for  numerical  computations 
than  A-11,  One  aethodj  therefore,  of  approximately  determin¬ 
ing  the  natriaant  is  to  break  up  the  matriaant  into  the  con¬ 
tributions  of  zones  as  in  A-12  and  then  treat  each  zone  as 
a  layer  of  constant  properties  so  that  A-14  can  be  used  for 
calculation.  This  is  essentially  the  layered  media  approach, 
but  the  use  of  the  matrlzant  calculations  is  especially 
handy  when  the  analytic  form  of  the  solutions  in  the  layers 
are  not  known.  The  convergence  of  the  layered  media 
approach  has  been  discussed  by  Ferraro  (Ferraro,  1962). 


As  we  mentioned  before^  a  waveguide  mode  problem  is 
usually  approached  as  a  characteristic  value  problem,  but 
in  the  Schumann  resonance  problem  we  are  dealing  with  essen¬ 
tially  only  one  waveguide  mode,  so  that  we  can  take  certain 
short  cuts.  If  the  one  mode  has  a  propagation  constant 


given  as 

k  =  cyf  -j- 

(per  unit  ©  ) 

^  vtv 

and  if  a  source  of  0 

is  imposed,  the  resulting 

waves  will  vary  as 

X  --  _e _ 


If  a  aolutlon  la  than  datarained  nuaarieallyj  the  phaaa 
of  tha  aoZutlon  allowa  ona  to  dataraina  tha  ratio  of 
and  tha  ratio  of  tha  anargy  fad  by  tha  aouroa  into  the  suida 
par  unit  d  to  tha  anargy  flow  along  the  guide  givaa  2  (3  . 


Appendix  B 


Two  PiMeneional  network  Solutlone 

The  si^pltet  difference  equetlon  epproxlaatlons  of  2nd 
order  partial  differential  equations  or  pairs  of  1st  order 
equations  lead  to  relationships  between  the  node  points  that 
can  be  interpreted  as  electric  network  equations.  A  great 
deal  of  literature  exists  on  nuaerlcal  isethods  of  solving  the 
equations  that  arise  from  potential  pdrobleas.  The  aost 
efficient  method  is  the  method  of  ovexrelaxation .  This 
method  is  still  very  efficient  when  dealing  with  dasiped 
electromagnetic  waves  in  conducting  media,  but  little 
analysis  exists  to  guide  one  towards  the  proper  amount 
of  overrelaxation  (  Cox,  I963  )•  Vhen  one  is  dealing  with 
essentially  propagating  waves  the  diagonal  dominance  of  the 
matrix  equations  is  lost  and  the  convergence  of  the  relaxa¬ 
tion  techniques  break  down,  and  one  must  consider  other 
techniques.  Vhen  one  looks  at  the  equations  as  network 
equations  one  is  naturally  led  into  a  technique  that  is 
straightforward  and  still  highly  efficient.  If  the  network 
under  consideration  was  a  one  dimensional  ladder,  no  electri¬ 
cal  engineer  would  ever  bother  with  simultaneous  solutions, 
but  instead  he  would  systematically  oospute  the  input  impe¬ 
dances  starting  at  one  end  and  working  back  to  the  beginning 
(jCron,  1959) •  The  partial  results  could  be  saved  and  from 
these  impedance  values  the  solution  to  any  given  source 


distribution  would  quickly  follow. 

The  vsry  ssas  stops  con  bo  tokon  for  two  dijssnsionol 
notworks^  but  now  ooeh  row  of  cAreuit  olonents  ploys  tho 
rolo  of  o  singlo  oloaont  in  tho  loddor  onolosy.  Itotrlx 
Invoroos  aust  thon  bo  perforaod^  but  tho  aotrix  sizos  ore 
deterained  by  tho  nusber  of  oloaonts  in  o  row  rothor  thon 
the  totol  nudber  of  eleaonts  in  tho  notwork. 

The  aethod  is  described  in  aotrix  longuogo  os  parti¬ 
tioning  and  triongularizotion  (Forsythe  and  Vosow«  I960). 

For  convenience  the  steps  in  tho  solution  olgoritha  oro 
repeated  below 

Given  a  network  and  current  soureoo^  X#  and  unknown 

voltages  V 

cr  =  Y  V  B-i 


which  we  partition  into  rows 


0  O 

o  -Ki  0 

O  ^  ”  ^Il3 


/  ^h2  *"  J 


b-4 


The  b's  are  the  adaittancee  of  the  circuit  eleaenta  connect¬ 
ing  the  node  points  to  the  nodes  of  the  rows  above  and 
below.  The  c's  are  the  adalttances  of  the  circuit  elements 
connecting  the  node  points  to  the  neighboring  nodes  on  the 
saae  row.  The  a's  represent  the  sum  of  all  the  admittances 
connected  to  the  nodcj  Including  the  admittance  from  the 
node  to  circuit  ground. 


Decompsaa  ^  Into  triangular  forms 

Y=  EF 


'too' 

'F,  G,  O 

f,  I  0 

0  I . 

f  - 

0  ^2 

/  '  " 

O  0  Fs 

% 

% 

• 

Sif 


E„F„=  C„ 

Gn  "  Gn 

En&rt 


B-7 


StartlJig  with  Fi  -  A|  th* 
bw  cA»tftln«d  by  th«  alsox^itlai 


•lawita  of  £  and  7  can 


Gn  = 


-I 


En  -  Cn  Ffi 

fhn  -  /\n-n”  En^n 


b-8 


LtX  FV  =  Z  = 


B-9 


Then  B-l  can  be  wrlttwi  aa 


El 


B-10 


Because  of  the  siaple  font  of  1  we  can  easily  obtain  Z  as 


B-11 


The  final  step  of  solving  for  \/  froa  ^  also  foUoirs 
directly 


-! 


The  last  eleaent 


The  rest  of  the  eleaants  follow  as 


B-12 


has  already  been  eoaputed  In  B-8 

The  eleaents  of  En  saved  as  a 

sort  of  Oreai's  function  for  solving  the  saae  network  with 
a  new  set  of  current  sources* 
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Appendix  C 

Dl»cwion  of  the  Sourcee  of  Srron  Iji  the  Cowputationa 

Ihere  are  three  type*  of  errors  Inherent  In  the  nuaeri- 
cal  results  presented  In  this  review.  The  first  type  re* 
suits  froa  errors  in  the  pliysical  laws  used  to  set  up  the 
theoretical  models .  The  only  significant  error  in  the 
theory  that  has  been  used  here  involves  the  simplification 
inherent  in  treating  particle  Interactions  in  teras  of  an 
effective  collision  frequency.  The  actial  particle  inter* 
actions  are  not  known  well  enough  to  allow  one  to  analyse 
the  effect  of  this  error  completely#  but  it  appears  that 
most  of  the  effect  can  be  taken  into  account  by  asking  the 
collision  parameters  slightly  frequency  dependant.  If  this 
is  the  case  the  effect  of  this  error  is  lost  in  the  second 
type  of  error#  that  due  to  improper  aodels  for  the  special 
distribution  of  the  important  physical  paraneters.  In  the 
long  run  the  errors  due  to  this  second  type  will  become  the 
data  which  will  be  used  to  correct  our  models.  This  cannot 
be  done#  however,  until  the  third  type  of  error  has  been 
eliminated,  that  due  to  the  aathematical  approximations 
used  to  obtain  numerical  results.  It  is  these  errors  which 
we  wish  to  re*examine  here.  These  errors  include  the  errors 
of  a  layered  approximation,  the  neglect  of  TE  modes  and 
higher  order  waveguide  inodes,  the  neglect  of  energy  return* 
ing  from  the  exosphere,  and  the  errors  of  the  network  approxi¬ 
mation. 
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It  18  known  that  th«  layered  media  approximation  conver- 
gea  to  the  correct  behavior  as  the  number  of  layers  becomes 
infinite^  so  that  an  empirical  method  of  testing  these  errors 
i6  to  carry  through  the  computations  for  thinner  and  thinner 
layers  and  stijdy  the  convergence.  In  our  computations  the 
layer  thicknesses  were  distributed  semi** logarithmically,  the 
layering  being  cruder  at  the  higher  levels.  Thus  a  nightime 
model  would  extend  to  300  KM  and  for  20  layers  the  lowest 
layers  were  about  5  KM  thick,  and  the  highest  layers,  30  KM 
thick.  The  ionosphere  above  300  KM  was  assumed  homogeneous. 
Each  layer  was  assumed  to  have  the  properties  associated 
with  the  middle  of  the  layer.  At  around  8  cycles  per  second 
even  a  20  layer  model  was  close  to  the  assymptotic  behavior, 
and  it  is  reasonable  to  expect  the  errors  are  within  about 
1%^  In  Table  C-1  the  behavior  of  the  solution  for  an  in¬ 
creasing  number  of  layers  is  shown. 


TABI£  C-1 


Effect  of  Number  of  Layers  on  Computed 
Propagation  Characteristics  Nightime  Model 


of  layers 

\X/q 

Q 

20 

.765 

2,44 

30 

.775 

2.28 

So 

.770 

2.32 

The  errors  involved  with  varying  the  height  to  the  last 
layer  were  of  a  comparable  magnitude.  In  the  daytime,  however, 
this  height  was  set  by  the  absorption  properties  of  the  iono¬ 
sphere,  If  the  height  was  extended  too  far,  the  solutions 
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beeaaa  unstable  baeausa  of  round«aff  errors  •  This  situation 
could  .be  aonitorad  from  the  partial  results. 

The  neglect  of  all  but  the  TBN  aode  was  a  basic  concept 
in  the  treataent  of  the  inhoaogeneoua  earth-ionosphere  wave¬ 
guide.  On  theoretical  grounds  this  was  well  justified  be¬ 
cause  of  the  severe  dai^ing  of  these  other  modes.  Some  check 
on  this  was  provided  by  the  cylindrical  models  which  did  not 
inhibit  these  modes.  Ihe  effect  of  the  higher  order  wave¬ 
guide  modes  could  not  be  discerned  in  the  solutions.  The 
presence  of  a  TE  mode  was  discernable*  but  its  energy  was 
also  insignificant. 

Since  30  per  cent  of  the  energy  escaped  out  of  the 
cavity  at  night,  and  since  the  aodel  did  not  account  for  any 
of  this  energy  returning,  some  concern  could  arise  about  the 
errors  of  neglecting  this  energy.  To  test  this  source  of 
error  models  were  set  up  where  this  energy  was  all  reflected 
back.  It  was  found  that  the  phase  velocities  were  changed  by 
as  much  as  4  per  cent  and  the  Q's  increased  from  3  to  30  per 
cent.  Since  only  a  fraction  of  the  energy  would  actually 
return,  and  since  the  coherency  of  this  energy  would  be 
lessened  by  variations  of  its  slow  path  through  the  exosphere, 
and  since  only  the  night  side  is  involved,  it  is  reasonable  to 
consider  this  error  as  not  more  lihan  1  per  cent. 

The  lumped  circuit  network  approximation  of  the  trans¬ 
mission  surface  equations  represented  a  siboificant  source 
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of  error.  For  a.  10*  x  20*  aodel  the  resonant  frequencies 
*fere  in  error  by  1,  2  and  3  psr  cent  for  the  first  three 
modes.  A  correction  factor  was  therefore  applied  to  the 
resonance  properties  that  were  coaputed.  ‘Die  aaplltude 
values  would  show  eve.i  larger  errors >  and  little  significance 
could  be  given  to  any  values  within  about  20  degrees  of  the 
source.  Again  correction  factors  could  be  devised  to  loprove 
the  results.  The  overall  confidence  In  the  muserlcal  re* 
suits  quoted  here  Is  about  coaperable  to  that  of  the  ezperl* 
mental  results  and  appears  to  be  adequate  for  the  present. 
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